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SUMMARY 


Centrifuge  experimental  techniques  provide  possibilities  for  laboratory 
simulation  of  ground  shock  and  cratering  effects  due  to  nuclear  weapons.  This 
premise  is  predicted  upon  the  results  of  a  similarity  analysis  which  indicates 
that  increased  gravity  is  a  necessary  condition  for  subscale  testing.  The 
objectives  of  this  investigation  were  to  examine  the  similarity  requirements  of 
this  type  of  subscale  testing  both  theoretically  and  experimentally.  To  do 
this,  a  series  of  centrifuge  experiments  were  performed  to  validate  theoretical 
similarity  requirements  as  well  as  to  determine  the  practicality  of  applying  the 
technique  to  dry  granular  soils  with  little  or  no  cohesion. 

Two  sets  of  experiments  were  performed.  The  first  was  a  series  of  ten 
shots  using  Ottawa  sand  as  a  convenient  and  wel 1-characteri zed  test  medium. 
Results  of  these  experiments  are: 

a)  Reproducibility  was  confirmed  in  the  centrifuge  environment. 

b)  Particle  size  effects  on  final  crater  configuration  were 
determined  to  be  negligible  for  the  cases  considered. 

c)  Validity  of  the  derived  similitude  requirements  was  demonstrated. 

d)  Scaling  rules  for  apparent  crater  size  and  equivalence  among 
different  explosive  types  were  derived. 

Upon  the  successful  completion  of  the  Ottawa  sand  experiments,  a  second 
set  of  experiments  was  undertaken  to  simulate  a  large-scale  cratering  event.  The 
JOHNIE  BOY  500-ton  nuclear  field  event  of  1062  was  chosen  for  the  centrifuge 
simulation  study.  The  prototype  geology  was  a  dry  homogeneous  alluvium  which 
was  straightforward  to  reconstitute  on  the  centrifuge.  It  was  a  wel  1-charac¬ 
teri  zed  event  for  which  considerable  computations  of  cratering  mechanics  were 
performed  and  it  was  used  as  a  standard  to  develop  a  full-scale  simulation 
method  for  cratering,  demonstrated  by  the  MINE  THROW  120-ton  Ammonium 
Nitrate/Fuel  Oi 1  (ANFO)  experiment  performed  in  1071. 


An  equivalent  full-scale  PETN  spherical  charge  configuration  was 
determined.  The  equivalQnce  criteria  was  based  upon  matching  the  kinetic  energy 


into  the  ground  as  well  as  the  shape  of  the  flow  field.  This  hypothetical  event 
was  then  simulated  at  subscale  on  a  centrifuge  at  345  G.  The  results 
satisfactorily  demonstrated  the  applicability  of  using  a  centrifuge  to  simulate 
a  small-yield  nuclear-cratering  event  in  alluvium. 

The  utility  of  the  centrifuge  method  is  based  upon  scaling  results  that 
indicate  the  equivalent  full-scale  explosive  yield  is  equal  to  the  actual  charge 
size  used  in  the  subscale  experiment  times  the  cube  of  the  centripetal 
accel  erati on. 
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SECTION  1 
INTRODUCTION 


1-1  BACKGROUND 

The  centrifuge  has  been  used  extensively  to  study  the  quasi-static 
response  of  soil  structures  to  lithostatic  loadings.  This  type  of  geotechnical 
modeling  has  gained  little  acceptance  in  the  United  States  but  has  been  widely 
used  in  Europe  and  in  the  Soviet  Union  during  the  past  thirty  years J 

The  applicability  of  this  method  to  dynamic  experiments  is  not  as  well, 

o 

if  at  all,  established.  Pokrovsky  and  Fyodorov-  discuss  experiments  designed  to 
study  the  effects  of  buried  explosives  used  in  the  construction  industry.  In 
chapter  eight  of  the  second  volime,  they  extensively  address  the  question  of 
ground  shock  as  applied  to  failure  prediction  of  soil  materials  for  excavation 
purposes.  In  volume  one,  a  few  pages  are  devoted  to  the  description  of  an 
investigation  of  the  soil  ejected  by  the  action  of  an  explosion.  These  early 
experiments,  perfumed  during  the  1940's,  used  1.5  gram  detonator  caps  as  the 
explosive  source  and  were  conducted  in  two  different  soil  media:  a  dry  sand  and 
a  moist  clay.  The  centifuge  was  operated  at  65  G  for  the  sand  and  at  111  3  for 
the  clay.  In  both  cases,  the  results  for  the  depth  and  the  radius  of  the 
apparent  crater  were  interpreted  as  satisfactory  when  compared  with  the 
calcul ated  full  size  dimensions  for  the  "natural  state  using  the  formula  of 
Boreskov ." 


Scott  and  Morgan1  s^  summary  remarks  concerning  the  entirety  of  the 
Russian  geotechnical  centrifuge  work  are  especially  appl ieeble  to  these  few 
reported  cratering  experiments.  He  states  that  "it  is  apparent  from  the  effort 
put  into  the  technique,  that  the  Russian  workers  consider  the  centrifuge 
technique  well  proven,  although  it  is  not  possible  to  discover  from  Pokrovsky's 
work  any  satisfactory  demonstration  of  the  correlation  between  model  and 
prototype  tests  for  any  of  the  studies  he  cites." 

The  only  other  reference  to  Soviet  cratering  experiments  in  an 

3 

accelerated  reference  frame  is  the  work  of  Viktorov  and  Stepenov.  Their  work 
addressed  soil  excavation  techniques  using  buried  explosives.  These 
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experiments,  conducted  with  charges  placed  at  optimum  and  greater  depth  of 
burial  in  moist  sand,  were  not  performed  in  a  centrifuge.  Since  they  ware 
modeling  the  effect  of  a  blast  on  the  throw-out  of  soil  and  rocks,  they  chose  to 
use  a  "linear  accelerator"  (presumably  a  rocket  sled)  to  eliminate  "the 
distortion  of  the  results  of  the  modeling  by  the  Coriolis  accel eration. . . . " 

In  a  recent  investigation,  Schmidt^  examined  the  application  of 
centrifuge  experimental  techniques  to  the  modeling  of  cratering  phenomena.  An 
oil  base  modeling  clay  was  used  to  investigate  depth  of  burial  effects  at 

centripetal  accelerations  up  to  430  G.  The  results  of  these  experiments  were  in 

3 

good  agreement  with  the  earlier  work  of  Viktorov  and  Stepenov,'  using  the 
results  of  a  dimensional  analysis  relating  soil  properties  and  explosive 
characteristics  for  the  different  experimental  conditions.  More  importantly, 
by  using  the  concept  of  a  gravity-scaled  charge  yield  parameter  derived  from  the 
dimensional  analysis,  the  results  of  both  sets  of  high-G  experiments  were  shown 
to  compare  favorably  with  data  obtained  from  the  Nevada  Test  Site  (NTS) 
cratering  series.  This  agreement  for  cratering  efficiency  as  a  function  of 
non-dimensional  depth  of  burst  demonstrated  the  relevance  of  using  subscale 
laboratory  test  performed  at  elevated  gravity  to  predict  large-scale  cratering 
events  in  the  field. 

This  work^  emphasized  cratering  efficiency  based  upon  apparent  crater 
vol  ime  but  left  observed  variations  in  crater  shape  unexplained.  In  addition, 
the  gravity  dependence  due  to  the  overburden  may  have  dominated  the  phenomenon 
for  buried  shots  which,  in  turn,  may  not  be  the  case  for  near  surface  events. 
Additionally,  the  use  of  a  cohesive  oil  base  clay  to  model  non-cohesive  soils 
may  have  introduced  offsetting  strength  effects  and  may  have  minimized  possible 
Coriolis  effects  due  to  the  cohesion  of  the  overturning  flap  during  crater 
formation.  These  concerns,  as  well  as  the  design  of  a  critical  experiment  to 
test  the  similitude  hypothesis,  provided  the  basis  for  the  program  described  in 
this  report. 

1-2  OBJECTIVES/ APPROACH 

The  overall  objective  of  this  program  was  to  demonstrate,  in  a 
quantitative  manner,  the  validity  of  using  a  1 abora tory- sc al e  centrifuge 
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experiment  to  simulate  an  actual  cratering  event.  The  JOHNIE  BOY  500 -ton 
nuclear  cratering  event  was  chosen  for  this  purpose.  It  was  conducted  in  the 
desert  alluvium  at  NTS  in  July,  1062. 

A  systematic  procedure  for  eliminating  experimental  uncertai nties  was 

devised.  A  new  rotor  for  the  Boeing  600-G  centrifuge  was  fabricated  to 

accommodate  non-cohesive  soils  utilizing  a  symmetric  swing-basket  design.  To 

verify  the  suitability  of  this  configuration,  ten  test  shots  were  perforned 

using  a  well-characterized  Ottawa  sand.  This  sand  tyoe  was  chosen  to  allow 

direct  comparison  with  the  large  body  of  high  quality  1-G  1 aboratory  cratering 

5  6 

data  generated  previously  by  Piekutowski.  ’ 

The  first  concern  was  reproducibil  ity  of  data  in  the  centrifuge 
environment  to  ensure  that  vibration,  windage,  slumping,  and  Coriolis  effects 
did  not  invalidate  proposed  crater  measurement  techniques.  The  next  step  was  to 
obtain  preliminary  confirmation  that  a  grav i ty- sc al ed  energy  concept  was 
applicable  to  a  surface  burst  configuration  in  Ottawa  sand.  This  will  be 
referred  to  as  the  constant  test  and  is  described  in  detail  in  Section  4-2. 
The  final  issue  addressed  by  this  series  of  tests  was  to  determine  the 
sensitivity  of  the  high-G  results  to  variations  in  sand- part icl  e-si ze  distri¬ 
bution.  The  results  of  these  Ottawa  sand  experiments  confirmed  the  suitability 
of  the  centrifuge  technique  for  this  type  of  modeling. 

As  is  shown  in  Section  2,  theoretical  requirements  for  non-trivial 
scaling  can  be  satisfied  with  a  model  material  in  the  centrifuge  identical  to 
that  of  the  prototype.  Therefore,  to  simulate  the  JOHNIE  ROY  event  in  a 
subscale  experiment,  it  was  imperative  to  use  a  soil  that  was  characteristic  of 
the  original  site.  A  suitable  material  was  obtained  from  the  Kirtland  APB 
(KAFB)  environs  and  supplied  by  R.  W.  Henny  of  AFWL.  Since  this  material 
differed  significantly  from  the  Ottawa  sand,  sample  preparation  techniques 
needed  to  be  developed  to  ensure  reproducibility.  The  constant  ^  test  was 
performed  on  samples  of  this  material  to  confirm  that  the  derived  similarity 
requirements  were  valid  for  the  KAFB  alluvium. 

For  chemical  explosives,  the  dependence  upon  different  explosive 
properties  can  be  accounted  for  by  including  the  charge  properties  in  the 
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dimensional  analysis.  For  nuclear  events  such  as  JOHNIF  BOY,  a  high 

explosive/nuclear  equivalence  must  be  determined.  To  accomplish  this,  an 

7  8 

analysis  was  performed  based  upon  the  so-called  MINE  THROW  technique.  ’  This 

involved  calculations  using  a  finite  difference  code  to  determine  the  size  of  an 

equivalent  full-scale  PETN  charge  size  and  an  associated  depth  of  burial.  The 

PETN  configuration  was  varied  until  the  resulting  early  time  flow  field  was  in 

g 

agreement  with  that  calculated  for  the  actual  JOHNIE  BOY  nuclear  event.  Tne 
results  of  these  calculations^  provided  a  hypothetical  full-scale, 

high-explosive  PETN  event  in  which  the  charge  radius  was  1.88  meters  and  the 
depth  of  burial  was  1.20  meters.  The  charge  mass  was  49.3  metric  tons  giving  a 
nuclear  equivalence  factor  of  13.4  percent.  These  calculated  results  were 

suppl  ied  by  41  len  J  * 

This  hypothetical  high-explosive  equivalent  event  was  to  be  simulated 
at  laboratory  subscale  using  a  1.25  gm  PETN  explosive  charge.  Based  upon  the 
similitude  requirement  g  E  =  constant  and  gL  =  constant,  the  experiment  was  to 
be  performed  at  345  0;  hence,  the  test  depth  of  burial  was  to  be  120  cm/345  or 
0.348  cm.  (Actual  placement  was  0.352  cm.)  To  determine  the  sensitivity  to  the 
calculated  equivalent  depth  of  burial,  Allen  suggested  that  a  second  shot  be 
performed  at  a  burial  depth  of  0.845  cm. 

A  total  of  six  shots  were  fired  in  this  simulation  series.  Two  shots 
under  identical  conditions  demonstrated  reproducibility.  The  next  two  shots 
using  two  different  charge  sizes  confirmed  that  the  KAFB  alluvium  satisfied  the 

constant  n?  test.  Two  final  shots,  performed  at  the  conditions  suggested  by 

Allen'*  ,  bracketed  the  actual  JOHNIE  BOY  results  as  predicted. 
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SECTION  2 

THEORY  OF  MODELING 


The  modeling  of  full-scale  cratering  events  by  subscale  centrifuge 
experiments  is  being  investigated.  To  validate  this  technique,  it  is  necessary 
to  show  under  what  conditions  a  subscale  experiment  does  simulate  a  full-scale 
prototype  and  to  determine  the  limitations  and  interpretations  of  subscale 
experiments.  To  examine  these  requirements,  the  relationships  between 
comparable,  but  different  scale,  experiments  are  derived  and  studied. 

Two  methods  are  commonly  used  to  derive  modeling  laws:  similarity 
analyses  and  dimensional  analyses.  Similarity  analyses  require  a  complete  set 
of  equations  adequate  for  describing  the  phenomena  in  question.  Dimensional 
analyses,  on  the  other  hand,  are  based  on  an  ad-hoc  choice  of  independent 
variables,  without  regard  to  the  physical  laws  relating  them. 

These  two  approaches  overlap  considerably.  The  governing  equations  of 
any  phenomena,  if  properly  posed,  must  be  expressible  in  a  dimensionally 
invariant  form.  Thus,  if  a  dimensional  analysis  is  performed,  using  those 
variables  occurring  in  the  complete  set  of  governing  equations  together  with 
appropriate  initial  and  boundary  conditions,  then  identical  results  are  obtained 
by  either  method.  Consequently,  the  identification  of  the  governing  equations 
serves  also  to  identify  the  controlling  variables. 

In  this  section,  a  similarity  analysis  of  cratering  phenomena  is 

1 9 

presented.  The  analysis  follows  along  the  lines  of  Crowley  c  and  Killian  and 
.  IT  . 

Germain  with  two  significant  differences.  The  general  balance  equations  of 
continuum  mechanics  are  employed,  as  opposed  to  particular  forms  used  in  the 
computer  codes  considered  by  Killian  and  Germain.  Secondly,  the  similarity 
requirements  of  the  general  balance  equations  appl icable  to  all  materials  are 
distinguished  from  the  requirements  imposed  by  particular  constitutive  equations 
used  to  describe  various  materials.  In  this  way  the  general ity  of  the  results 
is  more  apparent.  These  results  are  then  applied  to  the  problem  of  a  particular 
explosive  detonated  in  or  near  a  deformable  soil  medium  with  the  region  above 
the  soil  filled  with  a  gas  such  as  air.  It  is  assumed  that  each  of  the  three 
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media--the  soil,  the  explosive,  and  the  air,  can  be  modeled  as  a  continuum  with 
appropriate  constitutive  equations.  The  qeneral  field  equations  will  he 
discussed  first. 

2-1  SIMILARITY  ANALYSIS  OF  FIELD  EQUATIONS 

The  deformation  and  motion  of  any  continuum  must  satisfy  the  field 
equation  of  mechanics  in  their  qeneral  form.  The  qeneral  thennomechanica  1 

14 

response  of  materials,  qiven  by  Truesdell  and  Toupin,  include: 
hal ance  of  mass 

P  det  F  =  pQ  ( 1 ) 

balance  of  angular  momentum 

T=TT  (?) 

balance  of  linear  momentum 

div  T  +  pb  =  pa*  (3) 

and  balance  of  energy 

pe  =  Tr  (TL)  +  pr  -  div  q\  (4) 

These  equations  relate  the  following  seven  fundamental  fields 


x(X,t) 

position  vector 

(5.1) 

p(x\t) 

mass  density  per  unit  vol  ine 

(5.2) 

I(x,t) 

stress  tensor 

(5.3) 

b(X,t) 

body  force  vector  per  unit  mass 

(5.4) 

e(x\t) 

internal  energy  per  unit  mass 

(5.5) 

r(x\t) 

heat  supply  rate  per  unit  mass 

(5.5) 

q(X,t) 

heat  conduction  vector 

(5.7) 
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and  the  four  derived  fields 
v(X,t)  =  ^  x(X,t) 

a(X ,t)  =  ^  v(X,t) 

F(X,t)  =  Grad  x(X,t) 

t(x>>t)  =  [^t  [rV\t)] 


velocity  vector  (5.8) 

acceleration  vector  (5.9) 

deformation  gradient  tensor  (5.10) 

velocity  gradient  tensor  (5.11) 


all  expressible  as  functions  of  initial  position  X  and  time  t.  A  superposed 
arrow  denotes  a  vector  and  a  tilde  under  a  quantity  denotes  a  tensor.  The 
superscript  T  on  a  tensor  denotes  the  transpose,  and  the  inverse  is  denoted  by 
the  superscript,  -1.  The  operator  div  refers  to  the  divergence  with  respect  to 
the  spatial  position  x,  and  Grad  is  the  gradient  with  respect  to  X.  The 
operator  det  refers  to  the  determinate  and  a  superimposed  dot  denotes  the  total 
derivative  with  respect  to  time.  The  initial  mass  density  is  denoted  by  pq.  In 
addition,  these  fields  are  related  by  whichever  constitutive  equations  describe 
the  response  of  the  various  materials  and  by  the  appropriate  initial  conditions 
and  boundary  conditions.  Various  forms  of  the  constitutive  equations  are 
considered  later. 

At  surfaces  of  discontinuity  ( e . g . ,  shock  fronts),  these  equations  are 
augmented  by  jump  conditions.  For  a  shock  or  a  detonation  wave  with  local 
normal  n  moving  at  speed  U  into  an  undeformed  and  unstressed  material,  the 
following  equations  apply: 

ro^ss:  PQ  U  =  pU  -  p(^-n)  (6) 

momentum:  Tn  +  pQ  Uv  3  0  (/) 

energy:  PQ  U(e  -  eQ  +  ^  v  •  v  -  Q)  +  n  •  Tv  =  0  (8) 


16 


and  a  compatibility  condition 


U(F-I)n  +  v  =  0  (9) 

where  I  is  the  identity  tensor  and  the  variables  p,  v,  T,  e,  F  refer  to  the 
properties  behind  the  shock,  Q  is  the  energy  per  unit  mass  added  at  the  jump  for 
a  detonation  wave,  eQ  and  pq  are  the  energy  and  density  in  the  undeformed 
material  ahead  of  the  wave. 

Two  different  solutions  to  the  entire  set  of  eqs.  1  through  9  are  to  be 
compared.  All  of  the  quantities  associated  with  the  second  solution  will  be 
denoted  by  primes.  These  two  solutions  are  said  to  be  simi 1 ar  if  the  following 
rel ationshi ps 


x’(X'  ,t' )  =  «x  x()T,t) 

(10.1) 

P*(X\f)  =  ftp  p(X,t) 

(10.2) 

I* (x* ,t' )  =  nT  I(X,t) 

(10.3) 

b'(X' ,t')  -  ob  b(X,t) 

(10.4) 

e' (X* ,t’ )  -  «e  e(X,t) 

(10.5) 

r'  (X'  ,t' )  =  nr  r(X,t) 

(  10.6) 

q* (X*  ,f )  -  q(X,t) 

(10.7) 

Q'  (X’  ,t'  )  =  ntg  Q(X,t) 

(10.3) 

hold  at  homologous  points  defined  by 

x*  ■  <*x  x  (in. 9) 

for  homologous  time 

i 
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( in. in) 


t'  =  at  t. 

The  derived  fields  are  related  as  a  consequence  of  their  definitions  by 


v’(x\f)  *  ~  v(X,t) 
at 

(11.1) 

a'(X\t')  =  ^_a(X.t) 

<«t> 

(11.2) 

F'(X',f)  =  jjjj-  f(x,t) 

(11.3) 

L’(X’,f)  =  -L  L(x,t). 
t  - 

(11.4) 

Here  the  ten  various  quantities  a.  are  constants  called  scale  factors.  Note 
that  if  <*x  =  Qx  in  both  solutions,  the  reference  position  X  is  the  initial 
position:  x(t  =  0)  =  X.  This  is  henceforth  assumed,  leaving  nine  independent 
scale  factors.  Thus,  in  particular,  eq.  11.3  states  that  the  deformation 
gradient  must  be  the  same  at  homologous  points  and  therefore  the  strains  are 

identical  also.  In  more  detailed  terminology,  the  two  solutions  are  said  to  be 

15 

geometrically,  kinematically  and  dynamically  similar. 

The  question  of  the  existence  of  two  different  solutions  related  by 
these  similarity  requirements  is  to  be  investigated.  Both  sets  of  solutions 
must  satisfy  the  balance  equations  given  above,  regardless  of  the  constitutive 
equations.  This  requirement  will  give  certain  restrictions  relating  the  nine 
scale  factors.  Additional  restrictions  from  specific  constitutive  equations  are 
derived  separately  below. 

Assume  that  a  primed  solution  exists  and  substitute  the  primed  fields 
into  the  balance  equations  1  through  4  using  the  similarity  relationships  (eqs. 
10.1  -  10.10)  to  get 


r»  p  det  F  = 

p  i 


a  P„ 
p  o 
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(l?.l) 


F 


0 2.2) 


J  =  aj  I 


"t  .  -  «,  «x  + 

—  div  T  +  nt  a.  pb  =  -g- -  pa 

(x  -  p  b  2 

x  nt 


(X  nt 

P  e 


«•» 


pe  =  —  Tr(T  L)  +  m  a  pr 

„t  «t  -  °  r  ax 


div  q  . 


(13) 

(14) 


The  unprimed  fields  must  also  satisfy  the  balance  equations.  Hence  equations 
13.1  and  1 ? . 2  are  satisfied  identically. 


For  the  balance  of  linear  momentum,  eq.  3  can  be  used  for  div  T  in  eq. 
13,  giving 


a 


x 


aT 

- IX  IX. 

<*y  P  b 
.  J 


pb 


(15) 


4  . 

This  must  hold  for  all  X  and  t.  Thus,  unless  the  acceleration  field  a  is 

4  ^  4 

itself  a  scalar  multiple  of  the  body  force  b  or  either  a  or  b  is  identically 
zero,  it  is  necessary  that 


nt,  nt  nt 
-I  =  J2_  X 

"«  K)2 


and 


(16) 


a 


—  -  a  a. 
ax  P  b. 


(1/) 


The  most  common  body  force  per  unit  mass  is  a  constant  vector  proportional  to  a 
constant  gravity,  q.  In  this  case,  any  acceleration  field  which  is  a  scalar 

*  4 

multiple  of  b  is  at  most  a  rigid  body  motion.  If  b  ;  n,  then  eq.  17  is  not 

required.  If  a  ?  fl ,  as  in  statics,  then  eq.  16  is  not  required. 


Now  substitute  pe  from  the  balance  of  energy  eq.  4  into  eq.  14  to  get 


for  all  X  and  t.  Again,  assuming  the  various  fields  in  eq.  18  are  neither 
identically  zero  nor  simple  scalar  multiples  of  a  common  field,  it  is  necessary 
that 


a  (i  =  ol 
p  e 


(19) 


(20) 


a 

_£ 


a 


t 


(21) 


Requirements  brought  about  by  detonation  waves  in  the  explosive  can  now  be 
included.  For  this  case,  both  solutions  must  also  satisfy  the  jump  conditions 
across  the  wave  (eqs.  6-9).  An  analysis  identical  to  that  just  given  for  the 
balance  equations  produces  only  one  additional  restriction  which  involves  the 
scale  factor  <»q  for  the  specific  energy  of  the  explosive  Q 


a 


Q 


(22) 


Altogether  then,  the  balance  equations  plus  the  jump  conditions  provide 
six  similarity  restrictions  among  the  nine  scale  factors  relating  the  variable 
fields  as  defined  by  eqs.  10.1  through  10.10.  Eqs.  16  and  17  are  a  consequence 
of  the  balance  of  momentum,  and  eqs.  19,  20  and  21  result  from  the  balance  of 
energy,  and  the  energy  jump  condition  supplies  the  remaining  constraint,  eq.  22. 
The  balance  of  mass  and  the  balance  of  angular  momentum  provide  no  restrictions. 


and 


It  is  convenient  for  the  present  application  to  consider  ax,  at  and  «p 
as  independent.  The  six  restrictive  conditions  (eqs.  16-22)  can  then  be  used  to 
solve  for  the  remaining  scale  factors: 


a 


T 


% =  V(“t>‘ 


nt 

e 


"r  ’ 


(23.1) 

(23.2) 

(23.3) 

(23.4) 

(23.5) 

(23.6) 


Therefore  it  can  be  seen  that,  considering  only  the  balance  equations 
and  the  jump  conditions,  nontrivial  similar  solutions  are  allowed  with  arbitrary 
scaling  of  size,  time  and  density.  The  six  remaining  scale  factors  must  then 
satisfy  the  six  equations,  23.1  through  23.6. 

2-2  ANALYSIS  OF  SPECIFIC  CONSTITUTIVE  EQUATIONS 

In  addition  to  the  balance  equations  and  the  jump  equations  given 
above,  the  complete  solution  for  the  deformation  and  flow  of  any  continuous 
medium  depends  on  the  constitutive  equations  that  describe  the  behavior  of  that 
medium.  Therefore,  while  the  above  similarity  requirements  (eqs.  23.1  through 
23.6)  are  necessary ,  the  question  of  their  sufficiency  has  not  been  answered. 
Note,  however,  that  they  are  general  and  apply  to  all  materials. 
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The  nature  and  type  of  equations  that  describe  the  material  behavior 
for  a  given  medium  and  their  role  in  determining  similarity  requirements  is  now 
considered.  In  order  to  obtain  all  constraints  due  to  the  constitutive 
equations,  it  is  necessary  to  consider  various  types  of  such  equations  to  ensure 
that  the  complete  set  of  equations,  balance  and  constitutive,  give  a  well-posed 
problem  in  a  mathematical  sense.  That  is,  there  should  exist  unique  solutions 
when  appropriate  initial  and  boundary  conditions  are  given.  However,  the 
question  of  uniqueness  and  existence  of  solutions  to  this  complete  set  of 
equations  cannot  be  answered  at  any  level  of  generality.  Thus,  the  study  of 
well-known  classical  theories  is  utilized  where  experience  has  shown  that  unique 
sol utions  do  exist. 

A  first  special  case  considers  the  mechanical  deformation  of  a  linearly 
compressible  hydrostatic  medium.  All  thermodynamic  fields  are  omitted.  The 
stress  tensor  T  is  assumed  to  have  no  shear  stress  components  and  is  therefore 
described  by  the  pressure  p  and  the  identity  tensor  I, 

I  =  -p  I  (24) 

The  pressure  p  is  assumed  to  depend  linearly  on  the  density  change 

P  =  K0C|-P0/P)  (25.1) 

where  pq  is  the  initial  density  at  zero  pressure  and  KQ  is  the  bulk  modulus. 
Two  different  experiments,  designated  as  the  primed  and  the  unprimed,  are  to  be 
compared.  Each  of  these  may  be  in  a  different  material,  but  each  material  is 
modeled  by  the  linearly  compressible  equation  given  above.  Consequently,  for 
the  primed  material 

P'  •  Kq'  O-Pq'/p'  ).  (25*2) 

Are  these  two  constitutive  equations  (25.1  and  25.2)  consistent  with 
the  similarity  requirements  given  by  eqs.  23.1  through  23.6,  or  do  they 
introduce  additional  restrictions?  The  answer  is  easily  seen.  Since  the 
pressures  p  and  p'  are  related  by  the  parameter  which  satisfies  eq.  23.1,  it 


is  necessary  that  the  bulk  moduli  be  related  by  a-j-  also 


aT 


V  <1% V  .  „  /M2 

K0(l-p0/p)  K0  p  \at) 


(26.1) 


si  nee 


P0'  /P'  “ 


aP  po 

app 


P0/P 


(26.2) 


Therefore,  similarity  is  certainly  possible  in  this  case.  If  the  two  materials 
are  indeed  different,  and  the  two  experiments  have  different  distance,  tine,  and 
density  scale  factors,  then  the  bulk  moduli  of  the  two  materials  must  satisfy 
the  requirement  given  by  eq.  26.1.  For  two  experiments  conducted  in  the  same 
material  (a  =  1  and  Kq'  =  K  ) ,  eq.  26.1  can  be  satisfied  only  if  <*t  -  «x, 
leaving  only  one  independent  similarity  parameter,  the  size  scale  factor 
From  this  and  eqn,  23.1,  the  stress  and  the  strain  must  be  the  same  at 
homologous  points  in  the  two  experiments.  Hence  dynamically  similar  experiments 
of  different  size  scale,  arbitrary  ax  =  a  in  the  same  same  linearly 
compressible  hydrostatic  material  must  have  the  following  scale  factors 


r*T  =  1  (2/.1) 

ab=l/ax  (27.2) 

A  second  example,  a  calorically  perfect  heat-conducting  gas, 
illustrates  how  additional  constraints  are  introduced  based  upon  simple  thermo¬ 
dynamic  behavior.  The  stress  is  again  hydrostatic,  but  in  this  case  the 
pressure  is  given  by  a  perfect  gas  law 


P  =  (y-1)  oe 


(28) 


where  y  is  the  perfect  gas  constant.  The  heat  conduction  is  assumed  to  satisfy 
Fouri er' s  1  aw 


q  =  -k  grad  n 


(2d) 
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where  9  is  the  temperature  field  and  grad  denotes  the  gradient  with  respect  to 
the  spatial  coordinates.  Further  assume  that  the  internal  energy  e  and 
temperature  9  are  related  by  a  constant  specific  heat  C  as 

e  =  Cv  0.  (30) 

The  heat  conduction  can  then  be  rewritten  in  terms  of  the  conductivity  k  and  the 
gradient  of  the  internal  energy 


q  =  -  Jr-  grad  e.  ) 

Lv 

Equations  23  and  31  can  be  taken  as  the  fundamental  constitutive 
equations,  eliminating  the  need  to  consider  the  temperature  field.  Roth  eq.  28 
and  eq.  31  must  be  compatible  with  the  requirements  given  by  eqs.  23.1  through 
23.6.  Equation  23  together  with  eqs.  23.1  and  23.3  give 


(y  - 1 )  =  _P_  P'e 1  =  aPae  _  ,  (32) 

(V-l)  P'  p  e  aT 

Therefore,  similarity  is  impossible  unless  both  experiments  are  conducted  in 
materials  having  the  same  value  for  the  perfect  gas  constant-  This  is  a  well 
known  result  in  fluid  mechanics.^  Further  similarity  restrictions  on  material 
properties  can  be  obtained  using  eq.  10.7. 


n'  k'  y  1 
<1  -  -  7h  grad  e  = 

Lv 


aqq 


(33) 


where  grad'  denotes  the  gradient  with  respect  to  x' .  Using  eqs.  10.1,  10.5  and 
23.3 


grad' e'  ae  _ 

grad  e  '  «x  ^fc)2 

so  that  using  eqs.  31,  33  and  23.5 


(34) 


(35) 
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providing  an  additional  similarity  requirement. 


Equation  35  can  be  rewritten  as  follows 


(36) 


using  the  scale  factor  for  the  magnitude  of  the  characteristic  velocities  v  and 
v*  (eq.  11.1) 


and  the  scale  factor  for  characteristic  sizes  9  and  i' 


(37.1) 


Substituting  another  perfect  gas  relationship,  Cy  =  Cp/y  and  using  y  =  y‘ 
from  eq.  3?  gives 


pv£Cp  p'v'it'Cp' 

k  k1 


(38) 


This  ratio,  the  product  of  the  Reynolds  number  and  the  Prandtl  number  referred 
to  as  the  Peclet  number,^  must  be  equal  at  homologous  points  in  the  two  flows. 


The  sound  speeds  in  the  two  different  perfect  gas  media  are  given  by 


(39.1) 


and 


c’2  =  Y'  •  (39.2) 

Now  using  eq.  23.1,  the  similarity  requirement  for  stress 


and  the  scale  factor  for  density  ap,  the  ratio  of  sound  speeds  can  be  written 


XL  fil  £. 

Y  P  p' 


Using  eq.  36.1  for  the  velocity  ratio  gives 


(41) 


(42) 


Therefore,  eq.  32  which  requires  Y  =  Y1  leads  to  the  further  requirement  that 
Mach  number  be  equal  at  homologous  points.  In  this  example,  had  viscosity  been 
included,  separate  requirements  on  the  Reynolds  number  and  the  Prandtl  number 
would  have  also  been  obtained^1  in  place  of  the  requirements  on  Peclet  nimber 
given  by  eq.  38. 


A  further  reduction  occurs  if  it  is  now  assumed  that  each  of  the  two 
similar  experiments  is  in  fact  conducted  in  the  same  material.  Then  the  same 
material  properties  apply  for  each  material,  =  1 ,  and  eq.  36  for  material 

equivalence  gives 

“t =  K>2  (43) 

The  six  dependent  scale  factors  are  then  given  by 


(44.1) 

«b  =  1/(«x)3 

(44.2) 

%  -  1/(«X)Z 

(44.3) 

ar  ~  l/("x>4 

(44.4) 

,xq 

(44.5) 

nQ  =  i;(ftx)2 

(44.6) 

where  the  size  scale  factor  is  arbitrary.  Note  that  stress,  body  forces, 

internal  energy,  radiation  and  heat  conduction  all  scale  with  various  factors 
based  upon  the  size  scale  factor  a^. 


1 


i 

I 

A 

i 

j 

« 
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The  requirement  of  equation  44.6  implies  that,  in  experiments  with 
explosives,  the  heat  of  detonation  Q  cannot  be  the  same  in  two  similar 
experiments.  However,  consistent  with  the  assumptions  given  above,  and  insofar 
as  the  thermodynamical  constitutive  equations  stated  must  also  apply  to  the 
explosive,  the  two  experiments  must  have  the  same  thermodynamic  material 
properties  for  the  explosives.  These  two  requirements  are  in  practice  mutually 
exclusive.  This  conflict  arose  because  of  the  inclusion  of  heat  conduction 
effects.  That  is,  when  heat  conduction  effects  are  significant  in  experiments 
with  detonating  high  explosives,  similarity  between  experiments  with  different 
size  scales  cannot  be  achieved  with  the  same  materi al . 

Having  examined  these  two  special  examples,  more  general  considerations 
can  now  be  given.  Guided  by  the  last  example,  it  is  assumed  that  heat 

conduction  effects  can  be  ignored.  With  this  restriction,  consider  the  question 
of  the  complete  characteri zataion  of  a  medium  such  as  soil.  It  is  expected  that 
prediction  of  all  aspects  of  explosive  cratering  may  require  concepts  of 

compressibility,  nonlinearity,  yield,  fracture,  porosity,  cohesion  and  others. 

13 

As  stated  by  Truesdell  and  Noll  of  the  concept  of  dynamic  similarity:  "The 
more  complicated  the  constitutive  equation,  the  greater  the  number  of 
dimensionless  numbers  that  must  be  controlled  in  order  to  assure  dynamical 
similarity.  Ultimately  the  dimensionless  response  functional  itself  must  be  the 
same  for  the  two  materials  in  order  for  scaling  to  be  possible." 

In  agreement  with  this  statement,  a  third  important  special  case  is 

considered.  (Df  course,  as  the  previous  examples  using  simple  constitutive 
equation  show,  similarity  can  be  achieved  in  different  materials  if  the 

constitutive  equation  is  sufficiently  simple.  This  is  the  case  in  classical 
fluid  mechanics.  There  is  the  possibility  of  at  least  approximate  similarity  of 
some  aspects  between  different  real  materials,  even  when  their  total  behavior  is 
rather  complicated.)  Suppose  two  experiments  of  different  size  scales  are  to  be 
conducted  in  the  same  material  and  with  the  same  explosive.  In  this  case,  is 
similarity  possible?  The  answer  to  this  question  can  be  given  at  a  general 
level.  Consider  the  similarity  requirements  given  by  eqs.  23.1  through  23.6  but 
with  the  assumption  that  the  same  media  and  the  same  explosive  type  are  used  in 
both  experiments.  furthermore  assume  that  heat  conduction  effects  are  not 


significant.  Then  it  is  necessary  that  “p  =  1  »  “q  =  1  and  ^rom  eq.  23.6,  = 

o^.  Hence  there  remains  only  one  independent  scale  factor  giving 


«T  =  1 


a.  =  1/a 
b  x 


a  =  1 
e 


a  =  1  /a 
r  x 


aQ  =  1 


(45.1) 


(45.2) 


(45.3) 


(45.4) 


(45.5) 


Therefore,  similarity  is  indeed  possible,  as  long  as  the  constitutive 
equations  are  consistent  with  this  scaling.  The  stress,  the  strain,  the 
density,  the  internal  energy  and  the  heat  of  detonation  will  be  the  same  at 
homologous  points.  The  body  force  and  the  heat-rate-supply  term  must  scale  as 
the  reciprocal  of  the  size,  and  the  scale  factors  for  time  and  distance  arp 
equal.  Consequently,  if  «x  =  1/10  so  that  a  1/10  size  scale  experiment  is  to  be 
performed,  the  body  force  must  be  10  times  larger  as  must  be  the  heat-supply 
rate.  All  events  will  occur  in  1/10  the  time  over  1/10  the  distance,  and  all 
velocities  will  be  the  same. 

It  is  fairly  obvious  that  a  large  class  of  constitutive  equations  will 
be  compatible  with  similarity  at  this  level.  Any  relation  between  the  conserved 
quantities  such  as  the  stress,  strain,  and  internal  energy  will  be  directly 
compatible,  no  matter  how  complicated  or  nonlinear.  This  includes  nonlinear 
elasticity,  plasticity,  porosity,  spall  and  fracture,  the  Mi e-Grunei sen  equation 
of  state,  the  so-called  Jones-Wikins-Lee  (JWL)  equation  of  state  used  for  high 
explosives,  and  many  others.  Only  invariance  to  arbitrary  «  and  is 

necessary.  Thus,  it.  does  not  include  any  constitutive  behavior  not  invariant  to 
size  scaling  and  to  time  scaling.  For  example,  an  equation  between  stress  and 
rate  of  strain  would  not  be  consistent. 

This  observation  is  summarized  by  the  following  general  st  a*  "merit: 
Complete  and  exact  dynamic  similarity  can  he  achieved  between  two  Hi t *eren' 
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experiments  of  arbitrary  size  scale  in  the  same  material  as  long  as  the 
constitutive  behavior  is  1)  rate  independent  and  2)  has  no  inherent  size 
properties.  Of  these  two  requirements,  rate  independence  is  probably  the  most 
restrictive.  It  has  already  been  noted  that  this  rules  out  heat  conduction.  It 
is  interesting  that  the  property  of  rate- independence  emerges  as  the  significant 
property  that  allows  non-trivial  similarity  in  very  general  materials.  Whether 
this  restriction  is  of  practical  importance,  and  if  so,  for  which  materials,  can 
only  be  answered  by  experimentation. 

2-3  DIMENSIONLESS  PARAMETERS 

The  results  of  the  previous  sections  can  be  used  to  generate 
dimensionless  parameters  enabling  direct  definition  and  comparison  of  similar 
experiments.  In  particular,  the  tests  in  question  are  to  determine  the  final 
apparent  crater  formed  by  the  detonation  of  a  spherical  explosive  buried  in  a 
homogeneous  soil.  The  vol  ime  of  the  crater  is  determined  by  the  motion  history 
x(X,t).  This  history  is  itself  determined  by  the  complete  set  of  equations, 
balance,  jump  and  constitutive,  together  with  the  appropriate  initial  and 
boundary  conditions. 

It  is  recognized  that  the  behavior  of  the  soil  is  complex  and  requires 
complex  constitutive  equations.  Thus  it  is  probable  that  similarity  will  not  be 
achieved  unless  the  same  soil  is  used  for  similar  experiments.  With  this 
restriction,  as  shown  in  the  previous  sections,  similarity  is  possible  assuming 
only  that  the  constitutive  equations  describing  the  soil  are  independent  of  the 
scale  factors  for  size  and  for  time. 

A  very  qeneral  class  of  mechanical  constitutive  equations  are  those 
13 

which  Truesdell  and  Noll  have  called  simple  materials.  These  include  all 
materials  for  which  the  stress  tensor  at  any  material  point  at  the  present  time 
is  determined  by  the  past  history  of  the  strain  at  that  point.  This  includes 
all  types  of  non-linear  elastic,  el ast ic-pl astic ,  and  visco-elastic  materials, 
whether  solid  or  fluid. 

IB 

Truesdell  and  Noll  have  shown  that  the  complete  set  of  material 
constants  characteri zi ng  a  simple  material  can  involve  only  constants  that  are 
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dimensionless,  or  have  units  of  stress  or  time.  If  a  restriction  to  rate- 
independent  materials  is  made,  only  dimensionless  or  stress-unit  constants  can 
remain.  This  level  of  generality  still  includes  all  aspects  of  nonlinear 
elasticity,  rate  independent  elastic-plastic,  and  rate  independent  fracture  of 
any  type.  It  is  assumed  that  the  constitutive  equations  for  the  soil  are 
included  in  this  general  class  of  constitutive  equations. 

The  explosive  is  assumed  to  be  modeled  by  the  classical  Chapman-Jouqet 
theory,  where  the  combustion  products  behave  as  a  perfect  gas. 

With  these  assumptions  a  list  of  pertinent  parameters  can  be  given. 
The  explosive  behavior  is  determined  by 

Q  -  the  heat  of  detonation  per  unit  mass 
5  -  the  initial  density  of  the  explosive. 

The  soil  is  characteri zed  by 

p  -  the  initial  soil  density 
Y  -  a  material  strength  parameter. 

The  possibility  of  similarity  between  different  soils  will  be  included 
in  the  following  analysis.  As  discussed  above,  all  material  properties  of  the 
soil  are  either  dimensionless  or  have  units  of  stress.  The  material  strength  Y 
listed  here  is  assumed  to  have  stress  units.  The  inclusion  of  all  remaining 
material  constants,  whatever  their  nunber,  will  not  change  any  of  the  arguments 
to  follow  and  therefore  are  not  included.  The  perfect  gas  constant  y  for  the 
explosive  products  is  omitted  for  the  same  reason.  This  point  will  be  clarified 
sub sequently . 

The  ambient  air  is  not  considered,  and  the  initial  geometry  is 
determined  by 

a  -  the  explosive  charge  radius 
d  -  the  depth  of  burial 

In  addition,  the  solution  depends  on  the  body  force  arising  from 
g  -  gravity. 

The  scale  factors  for  all  variables  including  the  above  independent 
parameters  are  determined  by  the  three  independent  scale  factors  ,  a  and 
ap,  as  in  eqs.  23. 1-33.6.  These  three  scale  factors  can  be  expressed  in  terns 


of  the  actual  physical  parameters  of  the  two  experiments  to  be  compared.  For 
exampl  e 


a'  =  «x  a 

so  that 


Since  the  mass  W  of  explosive  is  given  by 


(46.1) 


(46. 2 ) 


W  = 


<Sa 


eq.  46.2  can  be  rewritten  as 


(46.3) 


(46.4) 


Likewise  the  mass  density  at  each  point  in  the  model  experiment  is  related  to 
the  full-scale  experiment  by  ap.  Hence  p'  =  apo  and  61  =  apS  f  giving 


« 


P 


(47.1) 


Hence  eq.  46.4  can 


(47.2) 


It  is  now  convenient  to  determine  <*  in  terms  of  the  specific  energy  O'  and  Q  of 
the  explosives  used  for  the  two  experiments.  In  particular, 


Q'  3  «Q  0 

so  that  using  eqs.  23.6  and  46.4 


(48) 


(49) 


30 


Solving  for  a 


t’ 


a 


t 


(50) 


The  scale  factors  for  the  other  variables  can  now  be  determined. 
Volimes  V‘  and  V  are  related  by  the  cube  of  the  scale  factor  for  distance  given 
by  eg.  47.2 


V*  =  (ax)2  V^V  (51) 

which  leads  to  the  definition  of  a  dimensionless  parameter,  referred  to  as  a 
"-group,  relating  the  variables  in  each  experiment 


V'P'  _  VP  _  (52) 

"W7  "  "W"  "  "1 


For  any  two  similar  experiments  "-|  must  have  the  same  value.  Scale  factors  for 
the  remaining  independent  parameters  of  interest  in  the  two  experiments  under 
comparison  can  also  be  expressed  in  terms  of  the  three  scale  factors  a^,  ap  and 
given  by  eqs .  47.2,  47.1  and  50.  Using  eqs.  10.4,  23.2  and  46.4,  the  scale 
factor  for  gravity  is 


g*  = 


V  = 


(53) 


This  leads  to  the  definition  of  a  second  "-group 


(54) 


The  depth  of  burial  for  the  two  experiments  satisfies  eq.  10.1 

vl/3 


axd  = 


W*G 

W  S' 


(55) 


so  let 


d* 


=  d 


(56) 


The  density  of  the  explosive  satisfies  eq.  10.2, 


31 


(57) 


6  =  a  6  =  6 

P  P 

as  used  in  eq.  47.1,  hence  let 


(53) 


The  material  strength  parameter  Y  has  stress  units.  As  a  consequence, 
it  must  transform  between  the  two  experiments  with  the  scale  factor  »-[■: 


Using  eqs.  46.4,  47.1  and  50 

r  6  ’  q  * 

Y  ”  6  Q 


(69) 


(60) 


The  appropriate  "-group  is  given  by 


V' 

S'Q1 


_Y_ 

<5Q 


(51) 


All  remaining  parameters  that  have  been  identified  are  either 
dimensionless  or  have  units  of  stress.  For  example,  dimensionless  constants 
include  the  perfect  gas  constant  y  of  the  explosive  products  or  dimensionless 
material  constants  for  the  soil.  These  dimensionless  constants  can  be  easily 
shown  to  be  required  to  have  the  same  value  in  similar  experiments.  Fach  such 
parameter  defines  an  additional  "-group.  Those  that  have  stress  units  will  lead 
to  "-groups  that  are  ratios  of  the  additional  parameters  to  the  variable  y 
identified  above.  Ultimately  the  discussion  will  be  restricted  to  experiments 
using  the  same  soil.  For  this  case,  all  additional  soil  material  property 
"-groups  are  constant  for  an  entire  series  of  experiments,  and  there  is  no  need 
t.o  include  them  specifically. 


The  rel ationships  among  the  scale  factors,  eqs.  ?3.1-?3.t,  applied  to 
the  eight  parameters,  seven  independent  and  one  dependent.,  of  a  cratering 
experiment  require  that  the  five  independent  "-groups,  through  ,  each  have 
the  same  value  for  similar  experiments.  Conversely,  similarity  between  two 
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experiments  will  be  achieved  when  these  five  "-groups  are  the  same  for  the  two 
experiments,  at  least  insofar  as  the  original  list  of  seven  independent 
variables  is  complete.  Any  additional  independent  variable  will  give  an 
additional  "-group  which  must  also  be  the  same.  These  "-groups  allow  the 
design  of  a  subscale  experiment  to  model  a  given  full-scale  event.  Arbitrary 
values  of  a ,  p  and  Q  can  be  chosen  or  equivalently,  values  for  the  scale 
factors  <*x,  ap ,  and  Then  there  is  a  subscale  similar  experiment,  having  the 

same  values  for  each  of  the  five  "-groups  as  does  the  full-scale  event. 

The  original  list  of  eight,  variables  is  not  independent.  The  vol  wie  V 
has  previously  been  identified  as  depending  upon  the  other  seven.  Thus  the 
following  relationship  can  be  written 

V  =  F(g,d,S,Y,r,a,Q) 

By  using  a  simple  change  of  variables,  eq.  f>2.1  can  be  rewritten  in  terms  of  the 
five  "-groups,  p,  a  and  t) 

it  i  =  ‘-y*  ( ^  2  ’ 71 3  ’  7I4  5  ’  *  ’  ^  ’  Q  ^ 

However,  as  noted  above,  a  given  full-scale  event  can  be  modeled  with  a  similar 
experiment  for  any  values  of  p,  a  and  Q  whatsoever.  As  a  consequence,  there 
exists  a  whole  family  of  similar  experiments,  where  p,  a  and  ')  may  take  on  any 
values  whatsoever,  and  for  which  through  are  the  same.  Thus  it  can  be 

concluded  that  p,  a  and  Q  can  be  varied  at  will,  but  as  long  as  through 
are  held  constant,  then  must  also  be  constant  and  depend,  at  most,  upon  the 
other  "-groups,  therefore  eg.  62.7  can  be  rewritten 

(^2  ,tt3 ,ti4  ,715  ^  (h2*3) 


where 


M.l) 

('>3 .?) 
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■3  ■  -  $r 

(63.3) 

tf4  =  p/6 

(63.4) 

Y 

n5  "  <S() 

(63.5) 

Therefore,  if  through 

"r  are  controlled  so  that  they  individually  have  the 

sane  value  in  two  experiments,  then  similarity  is  achieved  and  the  group  ^  will 
also  have  the  same  value.  In  this  way,  the  vol  ume  of  a  crater  in  a  large-  scale 
experiment  can  he  determined  by  performing  a  small-scale  similar  experiment,  and 

the  scale  factors  on  all  parameters  are  determined  by  the  choice  of  the  three 
parameters  0,  a,  0. 

^n  important  special  case  arises  when  experiments  are  conducted  in  the 
same  soil  and  with  the  same  explosive.  Then  both  and  are  constant  for  the 
experiments  as  would  any  additional  "-group  for  the  soil  and  the  explosive 
giving 


*1  =  W 

Having  picked  the  size  scale  factor  by  choosing  charge  size,  only  the  depth  of 
burial  d  and  the  gravity  g  need  to  be  controlled  in  order  to  ensure  similarity. 
For  surface  bursts  (zero  depth  of  burial),  "3  =  D  and 


This  relation  contains  all  the  information  about  the  volume  of  a  crater 
for  a  given  charge  in  a  given  medium.  Fixing  t.he  value  of  " uniquely 
determines  the  value  of  and  the  experiments  are  similar.  However,  to 

determine  the  function  KJ  ,  a  set  of  nonsimilar  experiments  wi'h  various  v  1 1  ies 
ot  must  he  performed.  It  is  important  to  note  that  ran  be  varied  either 
by  changing  the  charge  size  a*  Fixed  g  or  by  varying  g  for  a  fixed  charge 
rhat  is,  wifh  variations  of  gravity  alone,  using  a  rbirge  size  *0  sr’ 

laboritorv  experiments,  »he  dependence  of  vol  :rie  on  charge  s  ’ /e  (  <n  *>..  *  a  1  v 
de»  era  ’  ned  uS’ng  ‘he  ’-iroup  represent  a  *  i  on .  The  id  1  ’  ‘  '  'nil  |,> ;  e  n  t.  e  e- 
ran  1 1  Se  he  1e*  ec  "u>  1  *>v  • '  <  ;  e-- 1  nen*  S  of  ‘ais  t  yp*».  The  ■  ■  ■ . .  ...  ’  ,  ’  ..  . 
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to  other  dependences  upon  the  non-dimensional  groups  and  and  any  others  of 
conseouence.  This  illustrates  the  significance  of  being  able  to  vary  the 
gravity  field  strength  in  cratering  experiments. 


SECTION  3 

EXPERIMENTAL  TECHNIQUE 


3-1  CENTRIFUGE  DESCRIPTION 

The  Boeing  600  G  centrifuge  was  used  in  this  study.  This  machine  has  a 
dynamic  load  rating  of  60,000  G-kg  (66  G-tons)  at  620  rpm  and  was  constructed 
using  the  aerodynamic  housing  and  main  shaft  assembly  from  a  Gyrex  Model  2133 
centrifuge.  The  rotor  was  designed  and  fabricated  by  the  Boeing  Company  to 
incorporate  symmetric  swing  baskets  allowing  testing  of  non-cohesive  soil 
materials.  The  arm  radius  to  the  fully  extended  base  plate  is  139.7  cm  with  a 
maximum  payl oad  mass  of  260  kg  on  each  rotor  end.  An  overall  view  is  shown  in 
Fig.  1  and  the  details  of  the  swing  basket  and  soil  sample  container  are  shown 
in  Fig.  2. 

The  centrifuge  is  powered  by  a  30  horsepower  Eaton  Dynamatic  Model 
ACM-326-9 1 0B  drive  unit  incorporating  an  adjustable  speed,  constant  torque  eddy- 
current  clutch.  The  unit  also  has  electrical  dynamic  braking  allowing  shut-down 
from  maximum  rpm  in  less  than  30  seconds.  The  constant  output  speed  motor  and 
variable  drive  unit  are  shock  mounted  and  coupled  to  the  main  shaft  with  a  belt 
to  minimize  vibration. 

The  rotor  shaft  is  equipped  with  24  slip  rings  for  instrumentation 
channels,  three  220  V.a.c.  power  slip  rings  and  a  hydraulic  slip  ring  which  can 
accommodate  either  gas  or  liquid.  A  pair  of  motor  driven  Nikon  F2  35-mm  still 
cameras  are  hub  mounted  in  a  stereo  configuration.  These  cameras  provide  stereo¬ 
photo  coverage  of  one  rotor  end  with  a  maximum  framing  rate  of  six  per  second. 
Since  the  cameras  are  mounted  at  an  average  radius  of  13  cm,  they  experience 
centripetal  accel  erat.ions  of  up  to  60  G  at  maximum  rpm.  Semiautomatic  single 
framing  is  not  affected  by  this  loading,  however  in  the  continuous  motor  drive 
mode,  shutter  speeds  fall  out  of  calibration  above  350  rpn  (approximate  camera 
1 oad i ng  of  IB  G) . 

This  camera  installation  is  under  continuing  devel  opnent .  In  this 
program  it  provided  a  system  capable  of  dynamic  stereophoto  mapping  in  the  event. 


that  the  formed  craters  proved  to  he  unstable  due  to  slumping,  vibration  or 
windage.  For  the  soil  materials  tested  this  was  not  the  case  and  the  cameras 
were  used  to  confirm  earlier  tests  in  which  pre-formed  craters  were  spun  up 
successively  to  check  for  possible  shape  changes. 

A  data  reduction  technique  for  stereophotogrammetry  was  developed  for 
these  non-metric  cameras  under  static  conditions  using  object  space  calibration 
marks  on  the  soil  sample  container.  The  details  of  this  method  which  are 
directly  applicable  to  the  dynamic  system  are  given  in  Appendix  D.  Figs.  D1  and 
02  show  the  on-board  camera  installation  in  the  rotor  hub.  A  600  watt  quartz- 
halogen  lamp  provides  adequate  illumination  for  films  having  speed  of  ASA  125  or 
greater . 


3-2  EXPLOSIVE  CHARGE  DESIGN 

To  test  the  hypothesis  that  increased  gravity  could  be  used  to  scale 
energy,  various  explosive  charge  sizes  and  materials  were  required  for  experi¬ 
ments  to  be  performed  at  different  G  levels.  The  choice  of  explosive  device  was 
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influenced  by  the  work  of  Piekutowski ,  ’  who  reported  very  satisfactory  perfor¬ 
mance  from  precision  devices  supplied  by  the  R.  Stresau  Laboratories ,  Inc.,  of 
Spooner,  Wisconsin.  For  the  present  work,  four  different  charge  sizes  were 
employed.  These  included  the  two  basic  charge  sizes  used  by  Piekutowski 
(1.70-gm  lead  azide  and  0.49-gm  PETN)  and  two  larger  sizes. 

The  0.49-gm  composite  PETN  charge,  designated  CICS-5  by  Stresau 
consists  of  two  halves  of  a  concentric  sphere  of  PETN  pressed  around  a  silver 
azide  initiator,  designated  CISAS-5.  The  silver  azide  is  a  sphere  with  radius 
0.193  cm.  It  consists  of  a  maximum  amount,  of  0.130  grams  of  explosive  centrally 
initiated  using  a  notched  0.01?7-cm-diameter  tungsten  wire  to  form  a  spark  gap 
when  energized  by  an  electric  current.  The  actual  weight  of  PETN  in  this  device 
is  a  nominal  0.360  gram. 

The  two  larger  sizes  use  the  same  silver-azide  initiator  with  greater 
amounts  of  PETN  pressed  into  a  concentric  spherical  configuration.  The  larger 
of  the  two,  designated  CICS-4,  has  a  total  mass  of  4. OR  gm,  approx imately  an 
order  of  magnitude  greater  than  P iekutowsk i ’ s  basic  0.49-gm  PETN  charge.  This 
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device  contains  a  nominal  amount  of  3.960  gm  of  PETN  pressed  to  a  nominal 
density  of  1.70  gm/cc.  The  other  composite  charge,  designated  CICS-1.265 
contains  a  nominal  1.22  gm  of  PETN  and  has  a  total  mass  of  1.34  gm. 

The  actual  charge  weights  recorded  in  Appendices  A  and  B  are  those 
provided  by  the  R.  Stresau  Laboratory.  It  was  not  possible  to  confirm  these 
after  delivery  due  to  unknown  amounts  of  binder  and  surface  lacquer  that  are 

used  in  fabrication.  In  addition,  it  should  be  noted  that  the  weight  of  the 

silver  azide  in  the  detonator  could  vary  because  of  a  surface  condition  in  the 

mold.  The  silver-azide  weight  could  be  as  low  as  0.120  gm,  but  was  not  greater 

than  the  nominal  0.130  gm  specified.  For  the  purposes  of  calculations  it  was 
assumed  to  he  0.125  gm. 

Piekutowski*1  performed  a  series  of  calorimeter  tests  to  measure  the 
energy  of  detonation  in  order  to  arrive  at  a  charge  weight  of  PETN  which  would 
liberate  an  amount  of  energy  equal  to  that  of  the  1.70  gram  lead  azide  basic 
charge.  Bomb  calorimeter  measurements  performed  in  air  at  one  atmosphere 

provided  the  heat  of  combustion,  whereas  similar  tests  using  argon  provided  the 
heat  of  detonation.  His  data  obtained  for  the  heat  of  detonation  is  reproduced 
in  Table  1.  The  column  labeled  "Net  Heat  Release"  is  the  total  heat  release 
measured  using  the  argon  bomb  calorimeter  less  that  produced  by  the  silver  azide 
initiator  (calculated  to  be  2.34  x  10  ergs  based  upon  a  mean  weight  of  0.125  gm 
and  the  average  value  for  the  heat  of  detonation  of  silver  azide  given  by 

Piekutowski  to  be  1 .88  x  10  ^  ergs/gm).  Using  this  correction  for  the  composite 
PETN  charges,  the  average  heat  of  detonation  for  PETN  as  fabricated  in  this 
style  charge  is  5.66  x  1010  ergs/gm  with  a  coefficient  of  variation  of  7%. 

For  the  three  sizes  tested  by  Piekutowski  there  appears  to  be  some 
correlation  in  the  experimental  scatter  due  to  device  size.  This  is  to  be 

expected  since  there  is  a  certain  starting  distance  before  the  detonation  wave 
is  fully  established  in  the  PFTN.  However,  the  value  of  5.7U  x  10^  ergs/gm  at 
maximum  density  conditions  (1./7  gm/cc)  quoted  by  Piekutowski  is  only  2%  greater 
and  within  one  third  of  a  standard  deviation  of  the  value  calculated  from  his 
calorimeter  data.  The  mass  density  of  the  Stresau  charges,  reported  to  be 

approximately  1.7D  gm/cc  is  slightly  less  than  the  maximum  for  PI  TN  and  is 
consistent  with  the  value  calculated  above  of  5.66  x  10^'  ergs/gm. 
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The  scatter  in  the  energy  release  data  for  the  1.70  gram  lead  azide 
charges  has  a  coefficient  of  variation  of  1%  based  upon  the  four  air  shots  and 
two  argon  shots  in  the  bomb  calorimeter.^  The  value  of  7%  for  the  small  PETN 
charges  cited  above  is  probably  due  to  the  composite  design  employing  the  silver 
azide  initiator.  This  coefficient  of  variation  would  be  expected  to  decrease  as 
the  amount  of  PETN  in  the  composite  design  is  increased,  minimizing  the  effect 
of  a  finite  initiation  distance  in  the  PETN. 

These  calorimeter  results  can  be  used  to  calculate  the  effective 
density  6,  the  energy  per  unit  mass  Q,  and  the  energy  per  unit  volume  Qv  using 
the  following  definitions: 


X  _  3w 

(66) 

4iraJ 

Q  =  E/W 

(67) 

Qv  =  <$Q 

(68) 

where  W  is  the  total  charge  mass  including  the  initiator,  a  is  the  charge 
outside  radius,  and  E  is  the  total  heat  released  during  detonation  incuding  that 

due  to  the  initiator.  A  summary  of  these  material  properties  for  the  four 

charge  configurations  used  in  these  experiments  is  given  in  Table  2. 

For  surface  burst  experiments,  the  explosives  were  placed  at  zero  depth 
of  burial  by  carefully  excavating  a  void  equal  to  one-half  the  charge  volume 
using  a  micro-vacuum  consisting  of  a  piece  of  1 /B  i^h  diameter  Teflon  tubing 
taped  into  the  nozzle  of  a  standard  shop  vacuum.  The  electrical  leads  were 
taped  to  the  outside  walls  of  the  aluminim  soil  containers  and  strung  in  such  a 
way  that  no  interference  resulted  when  the  basket  swung  up  during  accel eration. 

The  charges  were  fired  using  a  standard  1 aboratory  d.c.  power  supply  which 

provides  up  to  40  A  at  40  V  through  the  electrical  slip  rings. 

The  charge  pi  acement  was  quite  stable  and  no  displacement  was  ever 
observed  prior  to  firing  when  on  several  occasions  the  centrifuge  was  shut  down 
from  full  speed  as  a  precautionary  measure  to  confirm  preshot  test  conditions. 
One  run  resulted  in  several  misfires  due  to  breakdown  in  the  enamel  insulation 
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where  it  ran  over  a  sharp  edge  that  had  formed  on  the  soil  container  due  to 
sample  preparation,  causing  the  charge  to  be  shunted  by  the  alumimm  structure. 
However,  the  stability  of  the  test  configuration  allowed  two  subsequent 
replacements  of  the  charge  on  the  "X"  end  of  the  rotor  and  one  replacement  on 
the  "0"  end  before  the  problem  was  uncovered.  In  all  cases  the  soil  showed  no 
evidence  of  being  disturbed,  except  for  minor  marks  on  the  surface  where  the 
lead  wires  touched  it. 

3-3  SAMPLE  PREPARATION 

Techniques  to  fabricate  the  Ottawa  sand  soil  samples  were  based  upon 
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the  prior  experience  of  Piekutowski ,  ’  who  visited  the  Boeing  facility  and 
participated  in  the  test  bed  preparation  for  shots  10-X  and  10-0  as  well  as  two 
prior  checkout  shots  using  detonator  caps. 

The  procedure  was  to  calibrate  the  volune  of  the  soil  container  using 
water  at  a  known  temperature.  This,  together  with  the  empty  weight  of  each 
container,  provided  a  basis  for  determining  the  average  density  of  the  finished 
soil  sample  nrior  to  placement  of  the  explosive.  The  sand  was  pluviated  by 
pouring  it  slowly  onto  a  sieve  and  allowing  it  to  free  fall  approximately  one 
meter.  The  chosen  sieve  size  would  just  allow  passage  of  the  larger  particles 
and  could  be  used  to  direct  the  placement  of  the  sand.  In  addition,  the  sieve 
was  continuously  shaken  from  side-to-side  to  further  disperse  the  sand  allowing 
it  to  fall  as  particles  instead  of  as  a  fluid  stream.  In  this  way  each  particle 
canes  to  an  equilibrium  resting  place  with  a  minimum  of  interference  due  to 
other  falling  particles.  This  technique  produced  the  maximum  obtainable  sand 
density  which  was  reproducible  for  a  given  sand  type  to  within  0.?b%. 

Three  different  sands  from  the  Ottawa  Silica  Company  of  Ottawa, 
Illinois,  were  used.  flintshot  was  desired  so  that  the  data  could  be  compared 
directly  to  that  of  Piekutowski.  Due  t.o  the  commercial  unavai  1  abi  1  i  ty  of 
flintshot,  the  first  two  shots  were  fired  in  Sawing  sand.  Later  a  supply  of 
flintshot  was  obtained  from  the  University  of  few  Mexico  Civil  [nginoerinq 
Research  facility  (CERf)  through  the  efforts  of  R.  W.  Henny  of  AfWL.  However, 
as  shown  in  Table  3,  the  size  distribution  of  Sawing  sand  is  quite  comparable  t.o 
that  of  flintshot,  albeit  a  bit  finer.  The  third  sand  used,  also  referred  to  in 


Table  3,  was  Banding  sand  which  is  considerably  finer  than  the  other  two  and 
provided  a  suitable  particle  size  variation. 

The  actual  measured  densities  obtained  for  the  various  Ottawa  sands 
were  about  1%  less  than  maximum  values  reported  by  Piekutowski  ’  in  his 

experiments  at  UORI.  This  probably  reflects  a  systematic  difference  in 
technique  between  the  two  facilities.  For  the  smaller  centrifuge  samples  total 
weight  could  be  obtained  directly  in  one  weighinq.  For  the  much  larger  liDRl 
samples,  multiple  weighings  were  required  and  any  spill-over  was  collected, 

weighed,  and  subtracted  from  the  total.  Also  a  mechanized  sand  conveyor  with  a 
rotating  squirrel  cage  was  used  to  disperse  the  sand  which  provided  very  uniform 
results.  Other  than  small  errors  in  the  actual  volume  determination  of  the  soil 

container,  the  most  significant  source  of  difference  is  attributable  to  the 

surface  leveling  technique.  For  the  centrifuge  studies,  a  piece  of  a  1  uni n in 
structural  angle  was  used  as  a  scree  to  smooth  and  level  the  surface.  This 
would  level  the  surface  to  approximately  +_  0.05  cm  which  is  the  size  of  the 
larger  sand  grains  which  build  up  on  the  metal  top  edge.  Generally,  two  or 
three  passes  with  the  scree  were  made  to  obtain  the  desired  surface  finish. 
This  is  in  contrast  to  a  single  pass  used  with  a  mechanized  scree  at  UORI. 

To  enhance  ejecta  definition,  the  surface  was  dusted  with  a  light  coat 
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of  black  spray  lacquer  as  suggested  by  Piekutowski.  Care  was  taken  to  avoid 
any  cementation  due  to  the  lacquer  application.  This  was  only  applied  to  the 
Ottawa  sand  samples  to  improve  contrast  on  the  otherwise  all  white  surface. 

The  alluvium  soil  posed  different  fabrication  problems  due  to  the  large 
percentage  of  fines.  The  above  technique  for  Ottawa  sand  placement  was  not 
applicable  at  all.  Instead,  soil  was  poured  into  the  container  in  3  to  4  cm 
thick  lifts.  After  placement  of  each  lift,  the  sample  container  was  dropped  a 
dozen  times  onto  a  platform  of  wooden  1 "  x  4  "s  from  a  height  of  approximately  ?'.) 
cm.  Greater  height  caused  the  soil  to  bounce  instead  of  compact. 

The  container  was  overfilled  t.o  a  height  of  approximately  3  cm  above 
the  top  edge.  At  this  point  both  samples  for  a  given  run  were  centrifuged  at.  <<:n 
run  (son  G)  for  in  minutes.  The  samples  were  then  removed  and  the  surface 
leveled  using  the  aluminiri  scree,  (or  the  nominal  four  percent,  moisture  content 


alluvium  (shots  1 S — X ,  13-0)  this  technique  produced  a  density  of  1.53  gm/cc. 
The  2. 65%  moisture  content  alluvium  (shots  17-X,  17-0)  measured  1.61  gm/cc. 

Samples  for  shots  16-X  and  16-0,  which  were  prepared  using  the  drop 
technique  but  not  centrifuged,  attained  a  considerably  lower  density  (1.45-1.40 
gm/cc)  which  proved  not  to  be  stable.  After  these  shots  were  fired,  the  ground 
zero  surface  had  subsided  approx imately  0.0  cm.  This  subsidence  probably 
occurred  in  two  stages,  some  during  spin-up  and  the  remaining  amount  during  the 
firing  of  the  explosive.  It  was  concluded  that  it  did  not  all  occur  during 
spin-up,  since  it  was  greater  near  the  crater  than  around  the  container  edge. 
As  a  result  these  data  points  for  apparent  crater  volume  were  ambiguous  due  to 
this  subsidence  of  the  ground  zero  surface  (the  maximum  lip  height  was  below  the 
original  level).  However,  they  satisfied  the  test  objective  to  ensure  the 
reproducibility  of  the  KAFB  alluvium.  Rased  upon  a  reference  plane  tangent  to 
the  minimum  surface  elevation  just  outside  the  crater  lip,  the  two  volumes 
differed  by  only  6". 

Sample  preparation  for  the  UDR l  1-G  control  shots  in  alluvium  differed 
in  that  the  centrifuge  could  not  be  used  to  obtain  maximum  density.  Instead, 
the  material  was  tamped  using  layers  of  approximately  3  cm  thickness.  In  this 
way  density  of  approximately  1.60  gm/cc  was  obtained.  This  procedure  would 
cause  a  tendency  toward  layering  of  the  sample,  but  according  to  Pi ekutowski ^ 
no  direct  evidence  was  observed,  although  the  scatter  was  much  greater  than  for 
his  previous  shots  in  Ottawa  sand. 

The  KAFB  alluvium  was  shipped  in  sealed  five  gallon  cans  lined  with 
plastic.  The  moisture  content  was  approx  imately  4°'.  for  all  the  cans.  Shots 
1 7 -X  and  17-0  were  fired  in  alluvium  that  was  exposed  to  laboratory  air 
coni  i t i tns  for  two  weeks.  This  material  had  Jr  ied  i  it  son ’  md  measured 
moi sture. 


All  the  soils  used  in  this  program  were  characterized  by  separate 
testing  performed  by  Shannon  and  Wilson,  Inc.,  of  Seattle,  Washington.  The 
results  of  their  tests  provided  failure  envelopes  for  confining  pressures  up  to 
?3.73  tons/ft  (37.56  bars)  and  the  associated  triaxial  loading  paths.  This 
data  is  given  in  Appendix  C.  These  material  properties  tests  provide  reference 
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conditions  for  the  soils  used  ind  this  i  nforri.it  ion  c.in  be  util'/ed  to  ■!■■  1  i" 
future  centrifuge  tests  on  these  sti>  aerials. 

1-4  ti:st  ^Ronmmi- 

The  symmetric  design  of  f  he  rotor  allows  firing  'wr.  shot  »,  Our  i  no  •  »„• 
same  run.  This  provides  an  experiment  al  convenience  hu'  is  no*  ne<  •*  ,sarv  *  or 

4 

dynamic  balance  as  originally  envisioned  since  neql  tuible  imnun*  s  : ! 

are  ejected  from  the  swing  baskets  'less  than  >c>  gm) .  dual  she's  •  h  • 

need  for  a  counterweight  and  do  provide  a  convenient  eon'ru]  in  crude 
preparation  as  well  as  test  conditions  and,  as  such,  were  ilwivs  employed.  :  .<  h 
shot  is  designated  by  the  run  niiober  followed  by  the  suf  ‘  i  <  '■  or  v  *  o  dt"i.v-  *  tie 
appropriate  rotor  end. 

The  entire  test,  sequence  following  t  he  mounting  ot  •  n>>  ♦  mis'e--  t , 
weighed  samples  is  sunnari  zed  by  the  following  procedural  <  bee*  list  ised  du'-mio 
the  run. 

(I)  Firing  circuit  continuity  check. 

(?)  Install  shorting  plugs  in  centrifuge  and  firing  umt(s). 

(3)  °lace  charges  in  soil  samples. 

(4)  Check  stray  voltage  in  firing  circuit. 

(5)  Connect  charge  leads. 

(b)  Close  centrifuge  aerodynamic  housing. 

(7)  Remove  shortinq  plugs  on  centrifuge. 

(?)  Close  safety  door  barrier. 

(0)  Remove  shorting  plugs  on  firing  unit.. 

(10)  Accelerate  rotor  to  fesired  rpn. 

(II)  Fi re  on  countdown. 

At  any  time  a  run  could  be  halted  and  restarted  a*  the  oiv  r  i  in  o' 
the  test  conductor.  This  occurred  on  many  oc:  r>  >ns  when  a  strty  wisher  ,i»- 
small  piece  of  wire  produced  an  explainable  sound  luring  s:un-up.  In  a  Id  •*  >on 
t.0  the  test  conductor,  there  is  a  centrifuge  operator  and  ,m  irlnime 
technician.  Each  operates  independently  and  is  responsible  for  r  he  si'"‘v 
his  aspect  of  the  procedure.  Upon  the  expiration  of  a  five  minute  pen-el  •  > 
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evacuate  *  hr*  e<  >'  isive  products  following  firing,  the  chamber  is  reopened. 
Stereo  ind  doc  orient  ary  photographs  i  '•»  t  tken  inned  lately  prior  to  any  handling 
of  •  he  seer  i  -tens  . 

All  th--  r im 1  >rs  produced  in  the  dry  soils  tested  in  this  program 

appeared  to  he  conpletel  v  ;  il>l «*  and  no  evidence  of  subsequent  s'npo  change  due 

•  •)  ;  ih  rat  ion,  windage  or  si  imp  inn  was  observed.  Therefore,  in  all  cases  the 

craters  co  i’ !  he  nr  *  wired  using  a  special  profiloi-fer  after  the  centrifuge  had 

stopped .  :’uns  !'i  •  h rough  1'  were  measured  using  the  depth  gauge  of  an  ordinary 

a  ■■  h  i  n  i  st 1 s  vernier  cal  i  per  in  coni  unr  t '  ;n  with  a  hori  zontal  reference  provided 
hv  i  -  -inch  s*et«l  mac  h  >  n  i  ■>*  1  s  scale  positioned  above  the  sample  using  gauge 

blocls.  These  "■*'  isurenen*  •»  were  ,ery  te<lious  and  gave  rise  K )  the  design  of  a 

pr  ‘  '  >”e‘er  uh  i  c  h  wts  f  tbrica*  *d  in  *ime  to  be  used  to  measure  the  craters  from 
runs  M  ‘hrough  I  .  Th‘>  ;.>■ )  f  i  1  onet  **r  is  a  scaled-down  ;.*'-si)n  or  one  originally 

L, 

des>gn.»  <  n  /  •  a*  >wsk  i  .  consists  of  Sr?  vent.  v-one ,  1  /,!.  inch  diameter  steel 

*•  »K  s ;  ■  r  ,i  ]  m  o,.  c  •  id  mig  r overage  of  4"’  cm  which  is  ■'P~  of  the 

\  ,  ■  ’  ”t‘n  1  (■"  •'*  iM'  . 

'be  r  i  i  ...»  !  wh  :  h  on*  u  *  s  *  he  soil  is  tapered  with  )  0 . 1 0  cm  flat  on 
i*.  'he  m  1  lengths  >u*  r  al  1  <•  1  so  t  ha*  1  i  s  pi  ac  >i  ten*  <  an  be  Measured  from 
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were  so  encourug i nj  that  additional  1-G  control  shit;  wore  made  with  the  larger 
PHTN  charges  designed  for  the  centrifuge  experiments  to  complete  this 
com pari  son . 

aiese  additional  shots  were  performed  by  A.  J.  Piekutowski  of  'JOT  I 

under  separate  contract  to  AFWL.  The  1-G  test  matrix  was  coordinated  by  R.  Vi. 

Henny  of  AFWL  and  the  explosive  charges  were  provided  from  the  lot  fabricated  by 

Stresau  Inc.  for  the  high-G  work.  The  details  and  results  of  these  experiments 

20 

were  supplied  by  Piekutowski  and  are  given  in  Appendix  B.  Tables  B3  and  B4 
give  detailed  data  for  the  larger  PFTN  charges  in  Ottawa  Flintshot  sand.  Table 
B5  provides  data  on  a  checkout  series  in  KAFB  alluvium  using  the  1.70  gram  lead 
azide  charge  design.  Tables  B6  and  B7  give  the  results  for  the  larger  »F.TN 
charges  in  KAFB  alluvium. 
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SECTION  4 

OTTAWA  SAND  RESULTS 


A  summary  of  the  ten  Ottawa  sand  centrifuge  shots  is  provided  by 

Table  4.  In  addition,  nine  1-G  control  shots  are  shown.  These  include  two  1.70- 

S  A 

gin  lead  azide  shots,  three  0.49-qu  PETN  shots,  and  four  shots  performed 

Of) 

expressly  to  evaluate  the  behavior  of  the  larger  PETN  charges  at  1  G. 
Complete  data  tables  as  well  as  stereo  photographs  and  crater  profiles  are  given 
in  Appendix  A  for  the  centrifuge  shots.  Appendix  B  contains  data  tables  for  the 
1 -G  control  shots.  All  of  these  shots  used  a  half-buried  spherical  explosive 
charge  configuration  in  a  homogeneous  soil.  For  all  of  these  zero-depth-of- 
burial  results,  the  crater  can  be  shown  to  depend  upon  only  a  single? 
non-dimensional  parameter,  referred  to  as  a  gravity-scaled  yield  parameter 

derived  in  Section  2.3,  quoted  in  eq.  54,  and  based  on  an  earlier  dimensional 

4 

analysis.  Application  to  these  experiments  is  discussed  below. 


4-1  PARTICLE  SIZE  DEPENDENCE 


Similarity  requirements  for  explosive  cratering  experiments  in 
identical  materials  are  satisfied  if  all  linear  dimensions  are  reduced  by  a 

factor  equal  to  the  reciprocal  of  the  gravity  field  strength  (eq.  45.2).  In  the 
case  of  a  granular  soil,  the  question  arises  as  to  whether  the  particle  size  is 
small  on  the  scale  of  the  experiment  allowing  the  soil  to  be  treated  as  a 

continuum  with  no  inherent  size  effect.  In  general,  since  the  particle  size 

distribution  determines  the  constitutive  behavior  of  a  given  soil,  a  conflict 
would  arise  if  particle  size  were  to  be  scaled.  This  is  a  non-trivial  question 
and  cannot  be  answered  in  general  for  a  material  without  appropriate  testing, 
since  the  possibility  of  an  inherent  size  property  (such  as  flaw  size,  or  pore 
size,  etc.)  can  control  the  phenomenon  under  investigation. 

Two  shots  were  conducted  to  pxplore  grain  size  effects  fo>" 

noncohesive  dry  sands.  Banding  sand  was  used  for  shots  P-X  and  13-y  which  wr» 
to  be  compared  with  shots  12-0  and  13-0,  fired  in  Hintshot.  Shots  12-0  and  12-X 
were  fired  at  451  G,  whereas  shots  13-0  and  13 -X  were  fired  at  3nn  <].  As  shown 
in  Table  3,  the  average  particle  size  of  the  Banding  sand  was  approximately  one- 
third  that  of  the  Tlintshot  sand.  1  arl  ier  wort  by  !'i  etutowsl  i  J  '  us  i  no  1  .  "> 


51 


✓ 


gram  lead  azide  charges  resulted  in  a  factor  of  2  difference  in  the  apparent 
crater  volume  obtained  for  the  two  materials.  For  these  1-G  tests,  the  crater 
vol  une  for  the  Banding  sand  at  maximum  density  (e  =  1.744  )  was  19?  cc, 
approximately  550|,»  less  than  the  FI  intshot  crater  vol  une.  The  corresponding 
crater  volumes  from  the  centrifuge  tests  for  both  the  306-G  shots  and  the  451 -G 
shots  were  only  6",  and  5%  smaller,  respectively,  for  the  Banding  sand  as  shown 
in  Table  5.  Note  that  the  maximum  density  obtainable  for  the  Banding  sand  in 
the  high-G  shots  was  somewhat  less  than  that  obtained  for  the  1-G  shot;  however 
the  former  was  a  stable  density  and  there  was  no  settling  observed.  This  lower 
density  could  reflect  a  slightly  different  size  distribution  since  the  UDR I  sand 
was  from  a  different  lot.  In  addition,  the  finer  sand  is  harder  to  pluviato 
since  it  is  less  uniformly  graded. 

An  observation  from  these  few  shots  is  that  for  increased  %  (large 
yield  or  large  gravity)  the  influence  of  particle  size  variation  on  crater  size 
decreases.  This  indicates  that  material  strength  effects  become  less  important 
for  large  Tr?.  Secondly,  it  confirms  that  the  Ottawa  sand  can  be  treated  as  a 
continuum  for  the  laboratory  scale  high-G  experiments.  This  also  suggests  that 
particle  grain  size  can  be  varied  to  achieve  sane  desired  behavior  without 
violating  the  scale  of  the  experiment.  For  the  case  of  saturated  sand,  perhaps 
grain  size  could  be  scaled  independently  from  experiment  size  to  obtain  the 
correct  pore  water  response. 

4-2  CONSTANT  TEST 

A  critical  test  has  been  defined  to  evaluate  the  applicability  of  using 
a  centrifuge  to  achieve  similarity  with  a  given  half  buried  explosive  in  a  given 
soil  material.  This  test  involves  measuring  crater  characteri sti cs  due  to 
charges  of  different  size  under  different  gravity  while  holding  the  value  of 
=  ( g/Q ) ( W/ « )  constant.  Then,  if  the  crater  characteri sti cs  for  a  soil 
material  are  functions  only  of  "j,,  the  material  models  used  in  deriving  the 
"-group  representation  are  adequate  (Section  2).  Furthermore,  this  allows 
charge  size  effects  to  be  determined  by  tests  in  which  gravity  alone  is  varied. 

For  "o  =  constant,  the  charge  weight  for  a  given  explosive  varies  as 
1/g  ,  fixing  the  relationship  between  W  and  g  for  fixed  "p .  To  test  this 
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relationship,  two  different  charge  sizes  of  the  same  type  explosive  charge  were 
fired  at  different  values  of  gravity  as  determined  by  the  ratio  of  charge 
masses.  When  the  test  matrix  was  devised,  it  was  assumed  that  for  the  same 
explosive,  charge  properties  would  not  depend  upon  size  and  therefore  the 
gravity  ratio  (g-j/g^)  for  a  given  charge  mass  ratio  (Wi/W^)  would  be  simply 


To  first  order,  this  assumption  is  valid;  but,  as  the  data  in  Table  2 
demonstrate,  the  presence  of  the  silver  azide  initiator  does  influence  the  value 
of  all  the  pertinent  charge  properties  for  the  composite  PTTN  charges.  Taking 
this  into  consideration,  the  prescribed  gravity  ratio  for  constant  " «  can  be 
expressed  as  follows 


Before  any  conclusions  regarding  the  application  of  centrifuge  results 
to  the  prediction  of  full-scale  1-G  events  can  be  made  for  a  given  material,  the 
constant  ^  test  must  be  shown  to  hold  for  different  charge  size  experiments 
performed  on  the  centrifuge  with  gravity  test  conditions  as  prescribed  by  eg. 
70. 

For  the  FI intshot  Ottawa  sand,  three  experiments  were  performed  to 
evaluate  compliance  with  the  constant  71  ^  test.  Shots  1 1  - X  and  11-0  were 
performed  at  451  G  to  confirm  reproducibility  using  the  nominal  1.34  gm  PTTN 
charges.  Shot  13-0,  a  4. OB  gm  PTTN  charge  fired  at  306  G,  provided  a  comparable 
value  of  for  which  the  cratering  efficiency  nv,  the  non-dimensional  crater 
radius,  11  ,  and  the  non-dimensional  crater  depth,  all  compare  favorably  with 
those  obtained  for  shots  1 1  - X  and  11-0.  The  slightly  larger  values  of  the 
dependent  variables  associated  with  shot.  13-0  are  consistent  wiith  the  lesser 
value  of  7 . 1 9  E  -  6 ,  as  opposed  to  7.54T-6,  for  the  shots  11-0  and  1  I  -  X . 

These  results  indicate  that,  to  well  within  experimental  scatter,  similarity  was 
achieved. 


4-3 


SCALING  RULES 


The  set  of  dimensionless  "-groups  provides  a  basis  for  comparing 
various  hi g h-G  and  1-G  results  at  differing  energy  levels.  As  the  previous 
analysis  shows,  experiments  at  a  constant  9ive  same  value  of  the  crater 
dependent  groups,  "  ,  "r  and  Thus  a  single  experiment  at  a  given 

furnishes  a  scaling  rule  for  all  similar  experiments.  However,  it  furnishes  no 
information  on  nonsiinilar  experiments,  such  as  varying  energy  and  hence  size  at 
fixed  gravity.  These  scaling  rules  are  discovered  experimentally  by  varying  "2 
in  a  series  of  experiments. 

Figures  3,  4,  and  3  show  each  of  the  dependent  groups  "  ,  "r  and 
respectively,  versus  "2  for  all  shots  included  in  Table  4.  As  can  be  seen,  all 
the  results  can  be  adequately  represented  by  a  power  law  fit,  such  as  the 
following  for  cratering  efficiency  "v: 

=  ky  =  const  (71) 

and  for  crater  radius, 

7TrTI2^  =  kr  =  const  (72) 

and  for  crater  depth, 

Y 

^h71?  =  kh  =  cons''  ( 73 ) 

These  then  are  the  form  of  the  function  r'J  as  given  in  (65). 


Expanding  eg.  71  using  the  definition  of  the  "-groups  gives  the 
explicit  form  for  crater  voltme: 


where,  from  the  experimental  results,  a  least  squares,  straight  line  fit  gives 
the  following  values  for  the  constants,  where  the  uncertainty  shown  for  the 
exponent  a  is  the  standard  error  of  estimate  (6 'IT  confidence): 


Surface  burst  scaling  rule  for  apparent  crater  volume  in  Ottawa  sand. 
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Figure  4.  Surface  burst  scaling  rule  for  apparent  crater  radius  in  Ottawa  sand. 


Figure  o.  Surface  burst  scaling  rule  for  apparent  crater  depth  in  Ottawa  sand 


V  =  (0.194  —  0.014)  (4 


1  / 3 -0.472  ±  0.005 


(74.2) 


Likewise,  from  the  best  fit  for  the  crater  radius. 


r  =  (0. 765  ±  0.030)  (4 


\l/3 


9./W 

Q  5 


1/3 


-0.159±0.003 


(75) 


and  the  best  fit  for  the  crater  depth. 


h  =  (0.1 54  ±  0.009) 


.1/3 


9. /W 

Q  6 


1/3  -i  -0.16410. 004 


(76) 


The  dependence  of  the  crater  dimensions  in  Ottawa  sand  on  all  of  the  independent 
variables  is  given  by  these  expressions. 

The  so-called  yield  exponents,  the  power  dependence  on  the  charge  mass 
with  all  other  variables  fixed,  are  determined  for  Ottawa  sand  from  these 
results  and  shown  with  uncertainties  corresponding  to  one  standard  error: 

(l-ci/3)  n  o no  +  n  on? 

V  «  W  =  W  ‘  u*uu^  (77) 


r  «  W0-P)/3  .  w0.280  ±  0.001 

h<*  Wn-Y)/3  .  w0.279  1  0.002 


(78) 

(79) 


If  9,  and  y  had  all  come  out  to  be  zero,  then  the  familiar  cube  root  scaling 
would  bp  applicable.  In  general,  this  is  not  the  case.  These  exponents  can  be 
considered  to  be  material  properties. 
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A1 ternatively,  consider  fixed  charge  size  and  type  at  varying  gravity. 
Then  eqs.  74,  75  and  76  reduce  to 


V  “ 


-a  -0.472  ±  0.005 

g  =  g 


(50) 


r  ct  g-P  =  g-0.159±  0.002 


(31) 


h«  g-T  =  g-0.164  ±  0.004  (32) 

giving  the  gravity  dependence  at  fixed  charge  size.  This  clearly  supports  the 
contention  of  the  previous  analyses  that  the  dependence  of  crater  formation  upon 
charge  size  can  be  determined  by  experiments  using  fixed  charge  size  but  varying 
g.  The  shot  10-X  at  463  G  using  1.70  gm  of  lead  azide  directly  simulates  a 
169-ton  mass  of  lead  azide  at  1  G  or  a  mass  of  2.86  kilotons  of  TNT  at  1  G. 

It  should  be  noted  that  the  variable  is  different  for  some  of  the 
tests  plotted  in  Figures  3,  4,  and  5.  However,  the  experimental  results 
indicate  that  Hy,  and  ^  do  not  depend  upon  this  parameter.  All  points  fall 

on  the  same  straight  line  when  plotted  versus  with  no  systematic  variation 
for  various  *4.  The  dependence  of  the  various  physical  parameters  actually 
varied  in  this  series  of  experiments  is  accounted  for  solely  by  those  occurring 
in 


In  summary,  eqs.  74,  75  and  76  are  referred  to  as  scaling  rules  for 
dynamically  nonsimi 1 ar  experiments.  They  are  functional  -el  at ionshi ps  for 
arbitrary  values  of  and  hence  for  arbitrary  size  experiments.  Such 

rel at i onshi ps  are  expected  to  have  different  functional  forms  for  different  soil 
materials.  These  scaling  rules  should  not  be  confused  with  scaling  relation¬ 
ships  based  upon  dynamically  similar  experiments  where  size  is  varied  while 
holding  the  value  of  as  well  as  other  pertinent  "-groups,  constant  .  In  the 
latter  case,  all  the  geometrical  features  are  identical,  differing  only  as  the 
size  scale  factor  In  the  former,  geometric  shape  variations  are  expected 

and  generally  are  the  case  as  " ?  is  varied. 
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The  important  observation  for  all  three  crater-size  parameters  shown  in 
Figs.  3,  4,  and  5  is  that  the  slope  of  the  line  determined  by  varing  charge  size 
at  constant  1-G  conditions  is  in  excellent  agreement  with  the  slope  of  the  line 
determined  by  varying  gravity  for  identical  charges.  This  provides  a  very 
convincing  argument  that  large  yield  can  be  simulated  with  small  charges  at 
elevated  gravity.  Furthermore,  it  confirms  the  use  of  n9  as  the  sole  parameter 
to  correlate  the  effects  of  gravity  and  charge  size  and  type  on  the  several 

forms  Try,  7ir  and  of  the  dependent  parameters.  The  parameter  ^  for  depth  of 

burial  was  held  constant  and  equal  to  zero.  No  data  were  obtained  on  the 
variation  with  respect  to  this  parameter. 

The  non-dimensional  form  of  the  gravity-scaled  yield  parameter  vorv 
adequately  accounts  for  the  large  differences  in  the  properties  of  the  various 
explosives.  For  the  composite  PETN  charges,  the  mass-averaged  values  given  by 
eqs.  66  and  67  were  used  for  specific  energy  Q  and  for  charge  density  6.  As  can 

be  seen  in  Table  2,  the  lead-azide  properties  and  the  PETN  properties  differ  by 

a  factor  of  4  on  specific  energy  and  a  factor  of  2  on  charge  density.  The 
various  values  for  the  four  different  explosive  charge  types  used  cover  the 
entire  range  encountered  for  common  chemical  high  explosives.  The  density 
parameter  was  varied  by  a  factor  of  2,  but  no  dependence  on  this  variation 
was  detected. 

This  result,  based  upon  the  utility  of  the  parameter,  provides  a 
basis  for  the  definition  of  an  equivalent  charge,  thereby  relaxing  the 
similarity  requirement  given  by  eq.  46.5 — that  the  same  explosive  type  is 
necessary  for  similar  experiments  in  the  same  soil.  This  has  a  particular 
significance  for  large-scale  simulation  using  the  centrifuge.  In  general,  the 
large  charges  used  in  the  field  consist  of  secondary,  explosives  which  are 
initiated  using  negligible  amounts  of  primary  explosives.  In  scaling  down  to 
laboratory  size,  the  amount  of  primary  necessary  for  initiation  becomes 
significant  in  relation  to  the  amount  of  secondary  changing  the  overall 
explosive  properties,  as  is  illustrated  in  Table  2.  On  the  basis  of  the  above 
results,  a  small-scale  equivalent  charge  of  a  different  type  can  be  designed  to 
simulate  the  large  prototype  explosive. 
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section  5 

J OHNIF  ROY  SIMULATION 


Six  high-G  shots  wore  performed  in  this  series.  The  first  pair  (ld-X, 
16-0)  confirmed  the  reproducibility  of  centrifuge  results  in  KAFB  alluvium. 
However,  in  spite  of  only  b"  variation  in  crater  volume,  considerable  subsidence 
of  the  ground  zero  surface  indicated  that  the  sample  fabrication  technique  was 
inadequate.  The  second  pair  of  shots  (17 -X »  17-0)  usinq  an  improved  sample 
preparation  method  as  discussed  in  Section  3-3  indicated  that  the  KAFB  alluvium 
satisfied  the  constant  "2  f  °st  of  Section  4-?  confirming  the  appl  icabi  1  i  t.y  of 
centrifuge  scaling  methods.  The  final  two  shots  (18-X,  18-0)  were  used  to 
provide  a  simulation  of  the  JOHN  IE  ROY  500-ton  event. 

Seven  1-G  control  shots  were  made  by  Piekutowsk  i  ^  to  provide 
additional  data  expanding  the  range  of  the  grav i tv- seal ed  yield  parameter, 

The  results  of  both  the  1-G  data  and  the  hiqh-G  centrifuge  data  are  stminarizod 
in  Table  6.  Detailed  records  of  the  six  high-G  shots  are  contained  in  Appendix 
A:  Tables  A7  through  A9  and  Figs.  A7  through  A9.  The  1-G  shot  records  are  in 

Appendix  R:  Tables  B5  through  B7  and  Figs.  R5  through  R7. 

5-1  KAFR  ALLUVIUM  SCALING  RULES 

Results  shown  in  Figs.  6,  7,  and  8  based  upon  t.he  limited  data 
generated  for  the  surface  burst  KAFR  alluvium  behavior  must  be  considered 
preliminary.  Of  the  six  high-G  shots  only  two  are  shown  ( 1 7 - X  and  1  / -0 ) .  The 
crater  shapes  determined  for  Ib-X  and  lb-0  were  based  upon  an  assumed  ground 
zero  plane  which  subsided  0.8  cm  during  firing.  Shots  1  ’ -X  and  18.-0  were 
performed  at  finite  depth  of  burial.  The  seven  1-G  shots  showed  considerable 
sender  indicating  possible  non-uniform  sample  preparation  as  discussed  in 
Section  3-3.  Further  tests  should  be  conducted  on  samples  prepared  in  'he 
centrifuge  and  fired  at  low  to  see  if  scatter  can  be  reduced. 

A  least  squares  straight-line  fit  »o  f he  data  indicates  that  *  he  slope 
is  considerably  less  *han  that  determined  for  Ot '  awa  sand  in  'he  previous 
section.  In  addition  tn  this  lesser  slope,  t.he  craters  are  smaller  '  han  their 

i 
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Figure  6.  Surface  burst  scaling  rule  for  apparent  crater  volume  in  KAFB  alluvium. 
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Figure  3.  Surface  burst  scaling  rule  for  apparent 


Ottawa  sand  counterparts .  For  ^  =  3.0F-8  (approximately  1-G  conditions)  there 
is  exactly  an  order  of  magnitude  difference  in  cratering  efficiency  between  the 
two  soils  as  shown  by  comparing  Figs.  3  and  6;  whereas  for  =  2.0E-5,  the 
cratering  efficiencies  are  equal.  Extrapolating  to  values  of  greater  than 
this  gives  larger  cratering  efficiency  in  the  KAFB  alluvium  than  in  the  Ottawa 
sand.  However,  these  observations  need  to  be  confirmed  by  additional  centrifuge 
shots  to  fill  in  the  range  of  *2  between  1.0E-7  and  1.0E-5  as  well  as  for  values 
greater  than  1 .OE-5. 

5-2  NUCLEAR-PETN  EQUIVALENCE 

It  was  shown  in  Section  2-2  that  centrifuge  experiments  rigorously 
satisfy  simil iarity  requirements  only  when  performed  using  identical  soil 
material  and  the  same  explosive  type.  As  a  practical  matter  the  above  test 
results  show  that  the  requirement  for  the  use  of  an  identical  explosive  can  be 
relaxed  to  the  use  of  an  equivalent  explosive  charge.  This  was  demonstrated  in 
Section  4-3  where  the  effects  due  to  different  charge  properties  were  accounted 
for  by  the  parameter. 

To  perform  a  centrifuge  simulation  of  the  JOHNIE  BOY  nuclear  event,  an 

equivalent  high  explosive  charge  was  determine^.  This  task  was  performed  by 

R.  T.  Allen  of  PacTech  under  separate  contract  to  the  Oefense  Nuclear  Agency. 

Allen  devised  a  hypothetical  full  scale  equivalent  PETN  charge  size  based  upon 

7  8 

the  so-called  MINE  THROW  technique.  ’  He  claimed  a  suitable  equivalence  based 
upon  two  criteria.  The  first  was  to  match  the  kinetic  energy  transferred  into 
the  ground  beyond  a  range  of  350  centimeters.  The  second  was  to  reproduce  the 
shape  of  the  flow  field  using  a  suitable  depth  of  burial  for  the  equivalent  PETN 
charge. 


The  energy  coupling  is  essentially  complete  after  the  first  few 

milliseconds  and  hence  can  be  conveniently  calculated  using  a  finite  difference 

wave  propagation  code.  The  standard  for  the  actual  JOHNIE  BOY  event  was  the 

9 

nuclear  calculation  performed  by  Maxwell  et  al.  Allen’s  results  indicated  that 
the  hypothetical  full  scale  charge  would  be  a  PETN  sphere  of  radius  138  cm 

buried  at  a  depth  of  120  cm.  This  corresponds  to  a  PETN  mass  of  49.3  metric 
tons  (49.3E5  gm)  with  a  total  energy  release  of  2.80E18  erg.  This  gives  a 
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nuclear  equivalence  factor  of  13.4%  based  upon  the  JOHNIE  BOY  nominal  energy 
release  of  500  metric  tons  of  TNT  nuclear  equivalent  (2.095E19  erg). 

The  calculated  energy  partitioning  for  the  JOHNIE  BOY  event  is  shown  in 
Fig.  9  along  with  Allen's  results  for  the  coupled  kinetic  energy  due  to  the  49.3 
ton  PETN  sphere.  The  best  equivalent  depth  of  burial  (120  cm)  based  on  these 
calculations  produced  a  flow  field  characteri zed  by  the  motion  of  the  reference 
hemispherical  surface  initially  located  at  a  radius  of  360  cm.  Its  motion  at 
3.6  msec  is  shown  in  Fig.  10  along  with  an  identical  reference  surface  from  the 
JOHNIE  BOY  calculation  for  comparison. 

Allen  chose  a  test  condition  of  345  G  based  upon  eq.  69  using  an 
explosive  charge  size  of  1.2  gm  PETN.  That  is 


Therefore  the  charge  was  to  be  buried  at 

(120  cm)  =  0.348  cm.  (84) 

A  second  shot  was  fired  on  the  same  centrifuge  run  as  a  backup.  For  this  it  was 
decided  to  increase  the  actual  depth  of  burial  to  approximately  1.5  times  the 
charge  radius  or  0.845  cm.  This  was  to  provide  a  determination  of  the 
sensitivity  of  flow-field  equivalence  to  depth  of  burial  and  to  bracket  the 
desired  result.  The  resuts  of  these  charge  equivalence  calculations  were 
presented  at  a  program  review  meeting^  prior  to  the  selection  of  the  actual 
alluvium  soil  material  which  was  eventually  used  for  the  simulation  experiments. 

5-3  EVALUATION  OF  RESULTS 

The  results  of  the  centrifuge  simulation  are  scaled  to  full  size  using 
the  similarity  relationship  given  by  eq .  45.2  as  g-jL-j  =  q^L-?.  These  art' 

compared  directly  to  the  actual  JOHNIE  BOY  crater  profile®  as  shown  in  Fig.  11. 
With  the  exception  of  the  aspect  ratio,  the  two  shots  18-0  and  18-X  bracket  the 
JOHNIE  BOY  crater  volume,  radius,  and  depth  as  anticipated.  The  aspect,  ratio 
was  slightly  larger  in  both  cases.  As  shown  in  Fig.  12,  the  results  ot  the 
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Figure  12.  Comparison  of  JOHNIE  bOY  crater  simulation  experiments, 
centrifuge  result  and  full  scale  MINE  THROW  event. 
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centrifuge  simulation  appear  to  be  equally  as  good  as  those  from  the  MINE  THROW 
event,'  a  full-scale,  high-explosive  field  event  which  used  120  tons  of  ANFO  to 
simulate  the  JOHNIE  ROY  crater.  It.  is  important  to  note  that  the  ANFO  charge 
was  not  spherical,  but  contoured  to  match  the  JOHNIE  BOY  pressure  history  at  a 
specific  isoharic  contour. 

The  centrifuge  simulation  of  a  shallow  buried  nuclear  event  was 
complicated  by  the  necessity  to  first  devise  a  chemical  explosive  equivalence. 
Therefore,  this  was  done  as  objectively  as  possible  by  an  independent 
calculation.  The  results  of  this  calculation  represented  the  design  of  a 
full-scale  simulation  of  JOHNIE  BOY  using  the  49.3-ton  sphere  of  nETN.  It  was 
this  hypothetical  event  that  was  simulated  at  small  scale  in  the  centrifuge. 
Uncertainties  in  the  nuclear  equivalence  could  be  responsible  for  any 
discrepancies  in  shape. 

Another  aspect  of  the  simulation  that  needs  to  be  considered  as  a 

possible  source  of  volume  discrepancy  is  the  requirement  for  reconstituting  a 

soil  sample  in  the  centrifuge  that  reproduces  the  average  _i_n  situ  properties  for 

the  JOHNIE  ROY  site.  JOHNIE  BOY  site  soil  density  was  reported  to  be  1 . 8S  gm/cc 

21 

based  upon  measurements  using  the  sand  displacement  method  in  a  test  pit.  The 

21 

moisture  content  at  the  2.5  foot  depth  was  3%.  Attempts  to  reproduce  this 
reported  value  for  density  with  the  KAFB  alluvium  test  sanples  were  not 
successful;  the  largest  attained  was  on  the  order  of  1.60  gm/cc,  a  value  more 
typical  of  desert  alluvium.  The  moisture  content  of  the  centrifuge  soil  sample 
was  4.1%. 

The  t_20%  difference  in  crater  volume  for  the  two  shots  is  well  within 

22 

the  combined  uncertainties  in  the  definition  of  the  JOHNIE  BOY  shot,  the  HE- 
equ i valence  deform inat ion,  and  the  experimental  accuracy. 

A  significant  shape  difference  between  the  actual  JOHNIE  ROY  crater  and 
the  centrifuge  craters  is  due  to  the  change  in  slope  halfway  up  the  rra*er  wall 
for  the  latter.  This  not.  only  leads  to  a  slightly  larger  radius,  bid  it  is  also 
noteworthy  thaf  the  lip  height  is  less  for  the  centrifuge  craters. 
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Three  factors  could  contribute  to  these  shape  differences.  A  spherical 
chemical  charge  might  not  adequately  model  the  late  time  low  pressure  expansion 
of  the  JOHNIE  BOY  event.  Another  possibility  is  that  the  model  soil  did  not 
possess  the  same  cohesion  near  the  surface  as  the  prototype  soil.  A  third 
possibil  ity  is  that  there  may  have  been  more  compaction  of  the  model  soil  near 
the  crater  attributable  to  lower  initial  density  than  that  reported  for  the 
prototype. 

It  is  interesting  to  note  that  the  surface  burst  shots  (centrifuge  runs 
16  and  17)  also  exhibited  the  discontinuity  in  slope  which  can  be  seen  quite 
prominently  in  the  photographs  given  in  Fig.  A7  through  A9.  The  UDRI  1-G 
control  shots  in  KAFB  alluvium  (also  surface  burst)  exhibited  a  shelf  at 
approximately  the  same  depth,  as  can  be  seen  in  Fig.  B5  through  B7,  and  the  lip 
heights  were  even  less  than  those  from  the  centrifuge  experiments. 

Considering  these  potential  sources  of  observed  discrepancies,  the 
agreement  of  the  two  centrifuge  shots  with  the  actual  JOHNIF  BOY  crater  as  shown 
by  Fig.  11  is  very  satisfactory.  These  results  strongly  support  the  theoretical 
premise  derived  in  Section  2  that  small-yield  shallow-buried  nuclear  cratering 
events  can  be  simulated  at  small  scale  in  the  laboratory  using  conventional 
explosives  in  a  centrifuge. 
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SECTION  6 
CONCLUSIONS 

(1)  A  similarity  analysis  was  performed  for  explosive  cratering  phenomena 

using  the  balance  equations  and  the  jump  equations  of  continuum  mechanics 
in  their  general  form.  Specific  constraints  on  centrifuge  modeling 

brought  about  by  various  types  of  constitutive  equations  representing 
material  behavior  were  identified  separately.  These  results  show  that  non¬ 
trivial  dynamic  similarity  can  be  achieved  using  a  centrifuge  with 
arbitrary  scaling  of  size,  density  and  time.  Furthermore,  complete  and 
exact  dynamic  similarity  can  be  obtained  between  two  different  experiments 
of  arbitrary  size  scale  in  the  same  material  as  long  as  the  constitutive 
behavior  is  rate  independent  and  has  no  inherent  size  properties. 

(2)  A  scaling  rule  for  surface  hurst's  in  dry  Ottawa  sand  was  developed  bast'd 

upon  high-G  experiments,  which  was  in  excellent  agreement  with  the  large 
body  of  1-G  work  performed  by  Piekutowski  .  The  crater  volume  depends 
upon  the  charge  mass  or  energy  (holding  all  other  variables  fixed)  to  the 
0.842  +^  0.002  power,  the  crater  radius  to  the  0.280  0.001  power,  and  the 

crater  depth  to  the  0.2/9  _+  0.00?  power  for  charge  mass  well  into  the 
kiloton  range. 

(3)  The  observation  for  all  three  crater-size  parameters  shown  in  Figs.  3,  4, 
and  5  is  that  the  slope  of  the  line  determined  by  varying  charge  size  at 
constant  1-G  conditions  is  in  excellent  agreement  with  the  slope  of  the 
line  determined  by  varying  gravity  for  identical  charges.  This  result, 
provides  convincing  evidence  that  larqe  yield  craters  can  be  simulated 
with  small  charges  at  elevated  qravity. 

(4)  A  non-dimensional  grav  i  ty- seal  ed  yield  parameter,  ij,  --=  g/0  (W/,1.)^^,  based 
upon  explosive  properties  and  qravity,  qave  quantitative  agreement  with 
four  different  composite  chemical  explosives  over  a  range  of  gravity 
between  1  G  and  451  G.  The  various  material  properties  for  the  difference 
charge  types  used  covers  the  on* ire  range  encountered  for  common  chemical 
high  explosives.  This  provides  a  basis  for  the  definition  of  an 
equivalent  charge.  It  further  suggests  the  possibility  that  a 
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nuclear-high-  explosive  equivalence  can  be  determined  directly  using 
centrifuge  experiments  in  conjunction  with  available  data  for  nuclear 
craters  • 


(5)  Preliminary  results  for  experiments  performed  using  KAFB  alluvium 

indicate  that  the  yield  dependence  is  considerably  different  than  that  for 
dry  Ottawa  sand.  This  scaling  behavior  needs  to  bo  character i zed  by 

additional  centrifuge  experiments. 

(6)  A  subscale  simulation  of  the  full-scale  500-ton  J0HN1E  BOy  crater  was 
satisfactorily  achieved  using  the  centrifuge  test  method.  A  nuclear-PETN 
charge  equivalence  was  calculated^  to  provide  a  hypothetical  full-scale 
equivalent  PFTN  spherical  center- initiated  charge  configuration  which  was 
then  simulated  at  subscale  on  the  centrifuge.  Two  centrifuge  shots 
bracketed  the  desired  result.  The  experimental  agreement  between  model 
and  prototype  for  the  size  of  the  apparent  crater  were  within  +17%  and 
-22%  on  volume,  >13%  and  -1%  on  radius,  +2%  and  -15%  on  depth,  and  +19% 
and  +11%  on  aspect  ratio. 

(7)  The  centrifuge  proved  to  be  an  effective  tool  for  modeling  explosive 
cratering  phenomena  in  dry  granular  soils.  No  measurable  Coriolis  effects 
were  observed.  Experimental  reproducibility  was  within  +_4%  on  volume  and 
+_2%  on  both  radius  and  depth.  Centrifuging  the  soil  samples  prior  to  test 
provides  an  effective  way  to  reconstitute  in  situ  material  properties  in  a 
reproducible  manner. 
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APPENDIX  A 

CENTRIFUGE  SHOT  RECORDS 


This  section  contains  all  the  data  for  the  elevated  gravity  centrifuge 
experiments.  A  table  is  given  for  each  run  which  includes  the  test  conditions 
for  the  two  shots  performed  at  the  opposing  rotor  ends.  Following  each  table  is 
a  figure  which  contains  a  comparison  plot  for  the  two  craters  as  well  as  a  set 
of  stereo  pairs  for  each  and  other  documentary  photographs.  Run  14  was  a 
demonstration  shot  using  modeling  clay  and  although  the  results  were  not 
discussed  in  this  report  it  was  included  for  future  referennce. 


PRECEDING  PAijE  BUNK -NOT  FliisED 
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Table  A1 .  Ottawa  sand  centrifuge  cratering  data. 


SHOT  NUMBER 

10-X 

10-0 

DATE 

9  Aug.  77 

9  Aug.  77 

PURPOSE 

charge  comparison 

charge  comparison 

CHARGE  DESCRIPTION 

C I LAS- 1 3  ( B-l ) 

CICS-5  (B-l) 

CHARGE  WT.  (gm) 

1.7000  PbN6 

.125  AgN3/ . 3601  PETN 

CHARGE  RADIUS  (cm) 

.508 

.391 

CHARGE  CONFIGURATION 

half-buried  sphere 

half-buried  sphere 

TEST  BED  MATERIAL 

Ottawa  Sawing  Sand 

Ottawa  Sawing  Sand 

TEST  BED  DENSITY  (gm/cc) 

1.774 

1.781 

MOISTURE  CONTENT 

nom  0  ' 

nom  0" 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

576 

576 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G's) 

463 

463 

CRATER  DIAMETER  (cm) 

7.74 

7.12 

CRATER  RADIUS  (cm) 

3.87 

3.56 

MAX  CRATER  DEPTH  (cm) 

.89 

.71 

CRATER  ASPECT  RATIO  (r/h) 

4.35 

5.01 

CRATER  VOLUME  (cc) 

27.05 

17.14 

LIP  DIAMETER  (cm) 

10.16 

8.64 

LIP  HEIGHT  (cm) 

.29 

.25 

LIP  VOLUME  (cc) 

-- 

MEAN  CRATER  PROFILE 

AVERAGE  OF  4  RADIALS 

FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

0.775 

0.673 

0.64 

0.883 

0.705 

1.27 

0.892 

0.711 

1.91 

0.845 

0.660 

2.54 

0.667 

0.470 

3.18 

0.451 

0.159 

3.56 

-- 

0.000 

3.81 

0.108 

-- 

3.87 

0.000 

-- 

4.32 

-- 

-0.245 

5.08 

-0.288 

-- 
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Figure  A1 .  (Continued,  crater  10-X). 


Figure  AT .  (Continued,  crater  10-0) 


Table  A2.  Ottawa  sand  centrifuge  cratering  data. 


SHOT  NUMBER 

11-X 

11-0 

DATE 

18  Aug.  77 

18  Aug.  77 

PURPOSE 

reproducibi 1 i  ty 

reproducibi 1 i ty 

CHARGE  DESCRIPTION 

CICS-1 . 265  (B-l ) 

CICS-1. 265  (B-2) 

CHARGE  WT.  (gm) 

.125  AgN3/l .2302  PETN 

.125  AgN3/l . 2248  PETN 

CHARGE  RADIUS  (cm) 

.565 

.565 

CHARGE  CONFIGURATION 

half-buried  sphere 

half-buried  sphere 

TEST  BED  MATERIAL 

Ottawa  Flintshot 

Ottawa  Flintshot 

TEST  BED  DENSITY  (gm/cc) 

1.776 

1.781 

MOISTURE  CONTENT 

nom  0% 

nom  OY' 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

568 

568 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G's) 

451 

451 

CRATER  DIAMETER  (cm) 

8.76 

8.62 

CRATER  RADIUS  (cm) 

4.38 

4.31 

MAX  CRATER  DEPTH  (cm) 

.99 

.97 

CRATER  ASPECT  RATIO  (r/h) 

4.42 

4.44 

CRATER  VOLUME  (cc) 

35.07 

37.95 

LIP  DIAMETER  (cm) 

11.56 

11.68 

LIP  HEIGHT  (cm) 

.33 

.33 

LIP  VOLUME  (cc) 

MEAN  CRATER  PROFILE 

AVERAGE  OF  4  RADIALS  FOR  EACH 

CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

0.991 

0.927 

0.64 

0.908 

0.965 

1.27 

0.775 

0.902 

1.91 

0.819 

0.927 

2.54 

0.794 

0.895 

3.18 

0.654 

0.705 

3.81 

0.394 

0.394 

4.38 

0.000 

-- 

4.57 

-- 

0.000 

5.78 

-0.330 

-- 

5.84 

-- 

-0.330 
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Figure  A2.  Comparison  of  crater  11-X  (above)  with  crater  11-0  (below) 
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Figure  A?.  (Continued,  crater  T  T -0 )  . 


Table  A3.  Ottawa  sand  centrifuge  cratering  data. 


SHOT  NUMBER 

1 2-X 

12-0 

DATE 

26  Aug.  77 

26  Aug.  77 

PURPOSE 

Particle  Size 

Particle  Size 

CHARGE  DESCRIPTION 

CICS-4  (B-2) 

CICS-4  (B-l ) 

CHARGE  WT.  (gm) 

.125  AgN3/3.9441  PETN 

.125  AgN 3/ 3 . 9600  PETN 

CHARGE  RADIUS  (cm) 

.826 

.826 

CHARGE  CONFIGURATION 

half-buried  sphere 

half-buried  sphere 

TEST  BED  MATERIAL 

Ottawa  Banding  Sand 

Ottawa  Flintshot 

TEST  BED  DENSITY  (gm/cc) 

1.677 

1.778 

MOISTURE  CONTENT 

nom  07. 

nom  0'' 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

568 

568 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G's) 

451 

451 

CRATER  DIAMETER  (cm) 

12.66 

12.82 

CRATER  RADIUS  (cm) 

6.33 

6.41 

MAX  CRATER  DEPTH  (cm) 

1.28 

1.31 

CRATER  ASPECT  RATIO  (r/h) 

4.95 

4.89 

CRATER  VOLUME  (cc) 

96.32 

101 .45 

LIP  DIAMETER  (cm) 

16.00 

16.00 

LIP  HEIGHT  (cm) 

.43 

.44 

LIP  VOLUME  (cc) 

65.18 

72.19 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS  FOP 

!  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  jcml 

0.0 

1.041 

1.283 

0.64 

1.191 

1.280 

1.27 

1.251 

1.283 

1.91 

1.172 

1.235 

2.54 

1.194 

1 .191 

3. 18 

1.114 

1 .159 

3.81 

1.041 

1.070 

4.45 

0.873 

0.918 

5.08 

0.641 

0.667 

5.  72 

0.327 

0.368 

6.35 

-i).  01  3 

0.  j25 

6.99 

-0.283 

-0.238 

7.62 

-0.425 

-0.406 

8.26 

-0.410 

-0.438 

8.89 

-0. 302 

-0.337 

9.  53 

-0.200 

-0. 260 

10.16 

-0.149 

-0.171 

10.80 

-0. 089 

-0.127 
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Figure  A3.  (Continued,  crater  l?-0). 
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Table  A4.  Ottawa  sand  centrifuge  cratering  data. 


SHOT  NUMBER 

1 3-X 

13-0 

DATE 

31  Aug.  77 

31  Aug.  77 

PURPOSE 

particle  size 

^2  fest 

CHARGE  DESCRIPTION 

CICS-4  (B-3) 

CICS-4  (B-4) 

CHARGE  WT.  (gm) 

.125  AgN3/3.9477  PETN 

.125  AgN3/3. 9491  PI 

CHARGE  RADIUS  (cm) 

.826 

.826 

CHARGE  CONFIGURATION 

hal f  buried  sphere 

half  buried  sphen 

TEST  BED  MATERIAL 

Ottawa  Banding  Sand 

Ottawa  Flintshot 

TEST  BED  DENSITY  (gm/cc) 

1.677 

1.782 

MOISTURE  CONTENT 

nom  0  , 

nom  O' 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

468 

468 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G's) 

306 

306 

CRATER  DIAMETER  (cm) 

13.30 

13.76 

CRATER  RADIUS  (cm) 

6.65 

6.88 

MAX  CRATER  DEPTH  (cm) 

1.56 

1.37 

CRATER  ASPECT  RATIO  (r/h) 

4.26 

5.02 

CRATER  VOLUME  (cc) 

119.30 

125.57 

LIP  DIAMETER  (cm) 

17.20 

16.40 

LIP  HEIGHT  (cm) 

.44 

.44 

LIP  VOLUME  (cc) 

65.47 

71.54 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS  FOR 

EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cni]_ 

0.0 

1.562 

0.978 

0.64 

1.467 

1.222 

1.27 

1.394 

1.340 

1.91 

1.359 

1  . 368 

2.54 

1.356 

1.330 

3.18 

1.299 

1.276 

3.81 

1.191 

1.235 

4.45 

1.038 

1.114 

5.08 

0.819 

0.876 

5.72 

0.498 

0.575 

6.35 

0.143 

0.260 

6.99 

-0.156 

-0.051 

7.6? 

-0.362 

-0.29? 

8.26 

-0.419 

-0.400 

8 . 89 

-0.359 

-0.416 

9. 53 

-0.254 

-0.346 

10.16 

-0.181 

-0.251 

10.80 

-0. 130 

-0.165 
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Table  A5. 


Permaplast"  clay  centrifuge  cratering  data. 


SHOT  NUMBER 

14-X 

14-0 

DATE 

19  Sept.  77 

19  Sept.  77 

PURPOSE 

tt 2  test 

tt 2  test 

CHARGE  DESCRIPTION 

C ICS-4  (B-5) 

CICS-1 .265  (B-3) 

CHARGE  WT.  (gm) 

.125  AgN3/3.9401  PETN 

.125  AgN3/l .2281  PETN 

CHARGE  RADIUS  (cm) 

.826 

.565 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

"Permaplast"  Clay 

"Permaplast"  Clay 

TEST  BED  DENSITY  (gm/cc) 

1.53 

1.53 

MOISTURE  CONTENT 

-- 

-- 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

471 

571 

GROUND  ZERO  RADIUS  (cm) 

124.5 

124.5 

CENTRIFUGAL  ACCELERATION  (G's) 

309 

454 

CRATER  DIAMETER  (cm) 

11.8 

8.34 

CRATER  RADIUS  (cm) 

5.91 

4.17 

CRATER  DEPTH  (cm) 

4.23 

2.59 

CRATER  ASPECT  RATIO  (r/h) 

1.40 

1.61 

CRATER  VOLUME  (cc) 

244. 

75.4  (+3) 

LIP  DIAMETER  (cm) 

14.4 

10.1 

LIP  HEIGHT  (cm) 

1.22 

.90 

LIP  VOLUME  (cc) 

79.3 

35.0 

MEAN  CRATER  PROFILE 

AVERAGE  OF  4  RADIALS 

FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

4.231 

2.594 

0.6 

3.986 

2.526 

1.2 

3.942 

2.397 

1.8 

3.826 

2.226 

2.4 

3.575 

1.923 

3.0 

3.196 

1.440 

3.6 

2.948 

0.774 

4.2 

2.366 

-0.045 

4.8 

1.590 

-0.749 

5.4 

0.800 

-0.824 

6.0 

-0.136 

-0.430 

6.6 

-0.829 

-0.046 

7.2 

-1.223 

-0.049 

7.8 

-0.859 

-0.059 

8.4 

-0.212 

-0.061 

9.0 

0.008 

-0.063 

9.6 

0.022 

-0.065 

10.2 

0.033 

-0.075 
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Figure  A5.  (Continued,  crater  14-X). 


Figure  A5.  (Continued,  crater  14-0). 


SHOT  NUMBER 

DATE 

PURPOSE 


Table  A6.  Ottawa  sand  centrifuge  cratering  data. 


1 5-X  15-0  i 

14  Oct.  77  14  Oct.  77 

seal i ng/reproducibil i ty  seal i ng/reproducibi 1 i ty 


CHARGE  DESCRIPTION 

CICS-5  (B-3) 

CICS-5  (B-4) 

CHARGE  WT.  (gm) 

.125  AgN3/. 3591  PETN 

.125  AgN3/.3598  PETN 

CHARGE  RADIUS  (cm) 

.390 

.390 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

Ottawa  Flintshot 

Ottawa  Flintshot 

TEST  BED  DENSITY  (gm/cc) 

1.784 

1.784 

MOISTURE  CONTENT 

nom  0% 

nom  0% 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

84.6 

84.6 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G’s) 

10 

10 

CRATER  DIAMETER  (cm) 

12.80 

12.80 

CRATER  RADIUS  (cm) 

6.40 

6.40 

MAX  CRATER  DEPTH  (cm) 

1.26 

1.27 

CRATER  ASPECT  RATIO  (r/h) 

5.08 

5.04 

CRATER  VOLUME  (cc) 

94.77 

95.18 

LIP  DIAMETER  (cm) 

15.60+ 

15.60+ 

LIP  HEIGHT  (cm) 

.35 

.33 

LIP  VOLUME  (cc) 

85.66 

78.87 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS 

FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

1.045 

1.130 

0.6 

1.089 

1.122 

1.2 

1.189 

1.206 

1.8 

1.239 

1.236 

2.4 

1.261 

1.267 

3.0 

1.221 

1.222 

3.6 

1.107 

1.107 

4.2 

0.933 

0.918 

4.8 

0.689 

0.686 

5.4 

0.412 

0.421 

6.0 

0.148 

0.147 

6.6 

-0.075 

-0.074 

7.2 

-0.259 

-0.252 

7.8 

-0.351 

-0.326 

8.4 

-0.348 

-0.320 

9.0 

-0.300 

-0.294 

9.6 

-0.232 

-0.217 

10.2 

-0.169 

-0.162 

DEPTH  (cm) 


Figure  A6.  Comparison  of  crater  15-X  (above)  with  crater  15-0  (below). 
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Figure  A6.  (Continued,  crater  15-X), 
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Figure  A6.  (Continued,  crater  1 5-0) . 
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Table  A7.  KAFB  alluvium  centrifuge  cratering  data. 


SHOT  NUMBER 

1 6-X 

16-0 

DATE 

3  Nov.  77 

3  Nov.  77 

PURPOSE 

reproducibi 1 i ty 

reproducibi 1 i ty 

CHARGE  DESCRIPTION 

C ICS- 1 .265  (B-9) 

C ICS- 1 . 265  (B-ll) 

CHARGE  WT.  (gm) 

.125  AgN3/l .2295  PETN 

.125  AgNo/l . 2353  PETN 

CHARGE  RADIUS  (cm) 

.565 

.565 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

KAFB  Desert  Alluvium 

KAFB  Desert  Alluvium 

TEST  BED  DENSITY  (gm/cc) 

1.449 

1.429 

MOISTURE  CONTENT 

4.1% 

4.1% 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

568 

568 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G's) 

451 

451 

CRATER  DIAMETER  (cm) 

9.40 

9.56 

CRATER  RADIUS  (cm) 

4.70 

4.78 

MAX  CRATER  DEPTH  (cm) 

1.91 

2.08 

CRATER  ASPECT  RATIO  (r/h) 

2.46 

2.30 

CRATER  VOLUME  (cc) 

48.11 

50.95 

LIP  DIAMETER  (cm) 

10.80 

10.80 

LIP  HEIGHT  (cm) 

.18 

.15 

LIP  VOLUME  (cc) 

MEAN  CRATER  PROFILE 

AVERAGE  OF  4  RADIALS  FOR 

EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

1.908 

2.075 

0.6 

1.942 

2.037 

1.2 

1.862 

1  .861 

1.8 

1.550 

1.512 

2.4 

0.898 

0.980 

3.0 

0.681 

0.669 

3.6 

0.434 

0.487 

4.2 

0.173 

0.223 

4.8 

-0.035 

0.008 

5.4 

-0.179 

-0.144 

6.0 

-0.142 

-0.108 

6.6 

-0.067 

-0.029 

7.2 

-0.027 

-0.008 

7.8 

-0.013 

0.000 

8.4 

-0.012 

0.006 

9.0 

-0.013 

0.000 

9.6 

-0.013 

0.003 

10.2 

-0.010 

-0.004 
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Figure  A7 .  (Continued,  crater 
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Table  A8.  KAFB  alluvium  centrifuge  cratering  data. 


SHOT  NUMBER 

17-X 

17-0 

DATE 

18  Nov.  77 

18  Nov.  77 

PURPOSE 

tt2  test 

n2  test 

CHARGE  DESCRIPTION 

CICS-1 . 265  ( B- 15) 

C ICS-4  (B-ll) 

CHARGE  WT.  (gm) 

.125  AgN3/l . 2281  PETN 

.125  AgN3/3. 9395  PETN 

CHARGE  RADIUS  (cm) 

.565 

.826 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

KAFB  A1 luvium 

KAFB  Alluvium 

TEST  BED  DENSITY  (gm/cc) 

1.607 

1.609 

MOISTURE  CONTENT 

2.65% 

2.65" 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

568 

468 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G's) 

451 

306 

CRATER  DIAMETER  (cm) 

8.44 

12.45 

CRATER  RADIUS  (cm) 

4.22 

6.23 

MAX  CRATER  DEPTH  (cm) 

1.35 

1.97 

CRATER  ASPECT  RATIO  (r/h) 

3.13 

3.16 

CRATER  VOLUME  (cc) 

30.23 

96.50 

LIP  DIAMETER  (cm) 

10.60 

15.00 

LIP  HEIGHT  (cm) 

.26 

.33 

LIP  VOLUME  (cc) 

21.34 

63.25 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS  FOR 

EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0 

1.227 

1.986 

.6 

1.355 

1.937 

1.2 

1.342 

1.947 

1.8 

0.917 

1.797 

2.4 

0.613 

1.547 

3.0 

0.488 

1.183 

3.6 

0.263 

0.954 

4.2 

0.013 

0.781 

4.8 

-0.201 

0.578 

5.4 

-0.250 

0.335 

6.0 

-0.161 

0.090 

6.6 

-0.094 

-0.139 

7.2 

-0.060 

-0.315 

7.8 

-0.042 

-0.339 

8.4 

-0.025 

-0.274 

9.0 

-0.016 

-0.209 

9.6 

-0.009 

-0.155 

10.2 

-0.008 

-0.113 

Figure  A8.  (Continued,  crater  17-0) 


Table  A9.  KAFB  alluvium  centrifuge  cratering  data. 


SHOT  NUMBER 

1 8-  X 

18-0 

DATE 

2  Dec.  77 

2  Dec.  77 

PURPOSE 

Johnie  Boy  Simulation  2 

Johnie  Boy  Simulation  1 

CHARGE  DESCRIPTION 

CICS-1 . 265  ( B- 17) 

CICS-1. 265  (B- 16) 

CHARGE  WT.  (gm) 

.125  AgN3/ 1 . 2243  PETN 

.125  AgN3/l . 2251  PETN 

CHARGE  RADIUS  (cm) 

.565 

.565 

CHARGE  CONFIGURATION 

.845  cm  DOB  (sphere) 

.362  cm  DOB  (sphere) 

TEST  BED  MATERIAL 

KAFB  Alluvium 

KAFB  Alluvium 

TEST  BED  DENSITY  (gm/cc) 

1.576 

1.570 

MOISTURE  CONTENT 

4.2% 

4.1% 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CENTRIFUGE  SPEED  (rpm) 

497 

497 

GROUND  ZERO  RADIUS  (cm) 

125 

125 

CENTRIFUGAL  ACCELERATION  (G's) 

345 

345 

CRATER  DIAMETER  (cm) 

12.16 

10.74 

CRATER  RADIUS  (cm) 

6.08 

5.37 

MAX  CRATER  DEPTH  'em) 

2.76 

2.26 

CRATER  ASPECT  RATIO  (r/h) 

2.20 

2.38 

CRATER  VOLUME  (cc) 

116.04 

77.74 

LIP  DIAMETER  (cm) 

15.00 

13.20 

LIP  HEIGHT  (cm) 

.37 

.32 

LIP  VOLUME  (cc) 

53.89 

35.18 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS  FOR 

EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

2.762 

2.257 

0.6 

2.760 

2.252 

1.2 

2.643 

2.201 

1.8 

2.352 

1.891 

2.4 

1.970 

1.454 

3.0 

1.449 

0.982 

3.6 

1.086 

0.761 

4.2 

0.852 

0.554 

4.8 

0.615 

0.312 

5.4 

0.333 

0.017 

6.0 

0.033 

-0.219 

6.6 

-0.218 

-0.322 

7.2 

-0.369 

-0.273 

7.8 

-0.365 

-0.171 

8.4 

-0.244 

-0.113 

9.0 

-0.162 

-0.074 

9.6 

-0.112 

-0.046 

10.2 

-0.069 

-0.023 
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Figure  A9.  Comparison  of  crater  1 8- X  (above)  with  crater  18  low  . 
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Figure  A9.  (Continued,  crater  18-X). 


Figure  A9.  (Continued,  crater  1B-0). 


APPENDIX  B 

UDRI  1-G  CONTROL  SHOT  RECORDS 


This  section  is  a  compilation  of  various  1-G  cratering  data  generated 

by  A.  J.  Piekutowski  of  UDRI.  The  data  for  the  larger  PETN  shots  in  both  the 

Ottawa  sand  and  the  KAFB  alluvium  was  generated  in  direct  support  of  this 

program  through  arrangements  with  R.  W.  Henny  of  AFWL.  The  data  for  the  1.7  gm 

lead  azide  and  the  0.4  gm  PETN  shots  in  Ottawa  sand  are  from  references  5  and  6. 

The  other  data  was  transmitted  as  a  personal  communication  from  A.  J. 
?n 

Piekutowski . 


PKECED1NG  PAOB  BLANK -NOT  FILMED 
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Table  Bl.  Ottawa  sand  1G  control  shots,  1.7  grams  lead  azide. 


SHOT  NUMBER 

DATE 

PURPOSE 

CHARGE  DESCRIPTION 
CHARGE  WT.  (gm) 

CHARGE  RADIUS  (cm) 

CHARGE  CONFIGURATION 

TEST  BED  MATERIAL 
TEST  BED  DENSITY  (gm/cc) 
MOISTURE  CONTENT 
TEST  BED  GEOMETRY 

CRATER  DIAMETER  (cm) 
CRATER  RADIUS  (cm) 

MAX.  CRATER  DEPTH  (cm) 
CRATER  ASPECT  RATIO  (r/h) 
CRATER  VOLUME  (cc) 

LIP  DIAMETER  (cm) 

LIP  HEIGHT  (cm) 

LIP  VOLUME  (cm3) 

MEAN  CRATER  PROFILE 


UDRI-9 
20  Jan.  72 
AFWL-TR-72-1 55 

C I  LAS- 1 3  (22) 
1.6961  PbN6 
.508 

half  buried  sphere 

Ottawa  Flintshot 
1.802 
nom.  0  a ; 
homogeneous 

20.4 

10.2 

2.29 

4.45 

436. 

260. 

.669 

321. 


UDRI-22 
9  Feb.  72 
AFWL-TR-72-1 55 

C I  LAS- 1 3  (106) 
1.7096  PbNg 
.508 

half  buried  sphere 

Ottawa  Flintshot 
1 .802 
nom.  0% 
homogeneous 

20.2 

10.1 

2.35 

4.30 

420. 

26. 

.763 

461. 


AVERAGE  OF  8  RADIALS  FOR  EACH  CRATER 


RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cnf 

0.0 

1.963 

1.950 

1.0 

1.988 

2.069 

2.0 

2.181 

2.306 

3.0 

2.288 

2.350 

4.0 

2.263 

2.300 

5.0 

2.131 

2.106 

6.0 

1.894 

1.825 

7.0 

1.513 

1.413 

8.0 

1.094 

0.975 

9.0 

0.550 

0.463 

10.0 

0.094 

0.038 

11.0 

-0.338 

-0.406 

12.0 

-0.581 

-0.650 

13.0 

-0.669 

-0.763 

14.0 

-0.575 

-0.675 

15.0 

-0.431 

-0.531 

16.0 

-0.338 

-0.431 

17.0 

-0.250 

-0.319 

18.0 

-0. 163 

-0.244 

19.0 

-0.088 

-0.181 

20.0 

-0.069 

-0.188 

21.0 

-0.050 

-0.175 

22.0 

-0.025 

-0.119 

23.0 

-0.013 

-0.106 

24.0 

0.000 

120 

-0.106 

Table  B2.  Ottawa  sand  16  control  shot,  0.4  gram  PETN. 


SHOT  NUMBER 

UDRI-146 

DATE 

24  Aug.  73 

PURPOSE 

AFWL-TR-74-182 

CHARGE  DESCRIPTION 

CICS-5  (103) 

CHARGE  WT.  (gm) 

.125  AgN ->/ .  3609  PETN 

CHARGE  RADIUS  (cm) 

.390 

CHARGE  CONFIGURATION 

half  buried  sphere 

TEST  BED  MATERIAL 

Ottawa  Flintshot 

TEST  BED  DENSITY  (gm/cc) 

1.796 

MOISTURE  CONTENT 

nom.  0 

TEST  BED  GEOMETRY 

homogeneous 

CRATER  DIAMETER  (cm) 

17.78 

CRATER  RADIUS  (cm) 

8.89 

MAX.  CRATER  DEPTH  (cm) 

2.02 

CRATER  ASPECT  RATIO  (r/h) 

4.40 

CRATER  VOLUME  (cc) 

284. 

LIP  DIAMETER  (cm) 

22.0 

LIP  HEIGHT  (cm) 

.581 

LIP  VOLUME  (cm3) 

216. 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS 

RANGE  (cm) 

DEPTH  (cm) 

0.0 

1.513 

1.0 

1.800 

2.0 

1.969 

3.0 

2.019 

4.0 

1.906 

5.0 

1.669 

6.0 

1 .319 

7.0 

0.888 

8.0 

0.406 

9.0 

-0.050 

10.0 

-0.394 

11.0 

-0.581 

12.0 

-0.544 

13.0 

-0.394 

14.0 

-0.231 

15.0 

-0.206 

16.0 

-0.119 

17.0 

-0.069 

18.0 

-0.050 

19.0 

-0.025 

20.0 

-0.019 

21.0 

-0.013 

22.0 

0.000 

Table  B2  (cont.)  Ottawa  sand  1G  control  shots,  0.4  gram  PETN. 


SHOT  NUMBER 

UDRI-11 7 

UDRI-1 20 

DATE 

11  June  73 

15  June  73 

PURPOSE 

AFWL-TR-74-182 

AFWL-TR-74-182 

CHARGE  DESCRIPTION 

CICS-5  (14) 

CICS-5  (21) 

CHARGE  WT.  (gm) 

.125  AgN3/.3602  PETN 

125  AgN3/ . 3602  PETN 

CHARGE  RADIUS  (cm>) 

.390 

.390 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

Ottawa  Flintshot 

Ottawa  Flintshot 

TEST  BED  DENSITY  (gm/cc) 

1.802 

1.796 

MOISTURE  CONTENT 

nom.  0 % 

nom.  0% 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CRATER  DIAMETER  (cm) 

17.5 

17.76 

CRATER  RADIUS  (cm) 

8.76 

8.88 

MAX.  CRATER  DEPTH  (cm) 

2.06 

2.13 

CRATER  ASPECT  RATIO  (r/h) 

4.25 

4.17 

CRATER  VOLUME  (cc) 

274. 

292. 

LIP  DIAMETER  (cm) 

22.0 

22.0 

LIP  HEIGHT  (cm) 

.569 

.581 

LIP  VOLUME  (cm3) 

217.0 

265. 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS 

i  FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

1.813 

1.925 

1.0 

1.938 

1.994 

2.0 

2.056 

2.125 

3.0 

2.038 

2.069 

4.0 

1.863 

1.950 

5.0 

1.625 

1.706 

6.0 

1.250 

1.313 

7.0 

0.825 

0.919 

8.0 

0.331 

0.406 

9.0 

-0.106 

-0.056 

10.0 

-0.431 

-0.425 

11.0 

-0.569 

-0.581 

12.0 

-0.513 

-0.556 

13. C 

-0.363 

-0.431 

14.0 

-0.263 

-0.288 

15.0 

-0.175 

-0.200 

16.0 

-0.125 

-0.156 

17.0 

-0.094 

-0.119 

18.0 

-0.044 

-0.081 

19.0 

-0.044 

-0.069 

20.0 

-0.031 

-0.031 

21.0 

-0.019 

-0.031 

22.0 

0.000 

-0.031 

23.0 

-- 

-0.031 

24.0 

-- 

-0.038 
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Figure  B2.  Comparison  of  craters  117,  120  and  146. 
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Table  B3.  Ottawa 

sand  1G  control  shots,  1 

265  grams  PETN.  ! 

SHOT  NUMBER 

UDRI-642 

UDRI-644 

DATE 

19  Oct.  77 

24  Oct.  77  ^ 

PURPOSE 

1G  Control 

1G  Control  1 

CHARGE  DESCRIPTION 

CICS-1 . 265  (B-5) 

CICS-1. 265  (B-4)  i 

*  CHARGE  WT.  (gm) 

.125  AgN3/1 . 2245  PETN 

.125  AgN3/l .2326  PETN  \ 

CHARGE  RADIUS  (cm) 

.565 

.565  j 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere  j 

s  TEST  BED  MATERIAL 

Ottawa  Flintshot 

1 

Ottawa  Flintshot  H 

TEST  BED  DENSITY  (gm/cc) 

1.802 

1.802  | 

MOISTURE  CONTENT 

nom.  O.Y 

nom.  0 l  P 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous  y 

CRATER  DIAMETER  (cm) 

23.6 

24.0  | 

CRATER  RADIUS  (cm) 

11.8 

12.0  <1 

MAX.  CRATER  DEPTH  (cm) 

2.50 

2.66  ij 

CRATER  ASPECT  RATIO  (r/h) 

4.72 

4.51  8 

CRATER  VOLUME  (cc) 

640. 

706.  1 

LIP  DIAMETER  (cm) 

30.0 

30.0  i 

LIP  HEIGHT  (cm) 

.838 

.863  1 

LIP  VOLUME  (cm3) 

798. 

814.  V 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS  FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

1.600 

1.650 

1.0 

1.681 

1.894 

2.0 

2.156 

2.319 

3.  C 

2.419 

2.569 

4.0 

2.500 

2.656 

5.0 

2.488 

2.613 

6.0 

2.319 

2.463 

'  |  7.0 

2.094 

2.206 

i  8.0 

1.731 

1.869 

9.0 

1.313 

1.481 

10.0 

0.838 

1.019 

11.0 

0.338 

0.488 

12.0 

-0.106 

0.019 

13.0 

-0.475 

-0.388 

14.0 

-0.719 

-0.681 

15.0 

-0.838 

-0.863 

16.0 

-0.763 

-0.844 

17.0 

-0.644 

-0.719 

18.0 

-0.531 

-0.581 

19.0 

-0.431 

-0.481 

,  20.0 

-0.363 

-0.381 

I  21.0 

-0.306 

-0.300 

1  22.0 

-0.263 

-0.256 

I  23.0 

-0.188 

-0.219 

1  24.0 

-0.156 

-0.169 
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Table  B4.  Ottawa  sand  1G  control  shots,  4.0  grams  PETN. 


SHOT  NUMBER 

UDRI-643 

UDRI-645 

DATE 

20  Oct.  77 

25  Oct.  77 

PURPOSE 

1G  Control 

1G  Control 

CHARGE  DESCRIPTION 

CICS-4  (BL-6) 

CICS-4  (BL-7) 

CHARGE  WT.  (gm) 

.125  AgN3/3. 9348  PETN  . 

125  AgN3/3. 9493  PETN 

CHARGE  RADIUS  (cm) 

.826 

.826 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

Ottawa  Flintshot 

Ottawa  Flintshot 

TEST  BED  DENSITY  (gm/cc) 

1.802 

1.802 

MOISTURE  CONTENT 

nom.  0"!. 

nom.  OT 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CRATER  DIAMETETER  (cm) 

32.6 

32.8 

CRATER  RADIUS  (cm) 

16.3 

16.4 

MAX.  CRATER  DEPTH  (cm) 

3.39 

3.54 

CRATER  ASPECT  RATIO  (r/h) 

4.81 

4.63 

CRATER  VOLUME  (cc) 

1677. 

1751. 

LIP  DIAMETER  (cm) 

42.0 

42.0 

LIP  HEIGHT  (cm) 

1.13 

1.17 

LIP  VOLUME  (cm3) 

1925. 

2062. 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADI ALS 

FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  XcmJ 

0.0 

2.400 

2.30C 

1.0 

2.338 

2.388 

2.0 

2.638 

2.869 

3.0 

2.931 

3.194 

4.0 

3.156 

3.406 

5.0 

3.331 

3.513 

6.0 

3.388 

3.538 

7.0 

3.356 

3.488 

8.0 

3.244 

3.294 

9.0 

3.031 

3.069 

10.0 

2.750 

2.794 

11.0 

2.394 

2.475 

12.0 

2.019 

2.125 

13.0 

1.556 

1.656 

14.0 

1.081 

1.163 

15.0 

0.600 

0.638 

16.0 

0.119 

0.163 

17.0 

-0.281 

-0.256 

18.0 

-0.638 

-0.625 

19.0 

-0.919 

-0.925 

20.0 

-1.081 

-1.106 

21.0 

-1.125 

-1.169 

22.0 

-1.050 

-1.119 

23.  C 

-0.950 

-1.000 

24.0 

-0.800 

-0.856 
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Table  B5.  KAFB  alluvium 

1G  control  shots,  1, 

.  7  grams  lead  azide. 

SHOT  NUMBER 

UDRI-646 

UDRI-647 

DATE 

11  Nov.  77 

14  Nov.  77 

PURPOSE 

1G  Control 

1G  Control 

CHARGE  DESCRIPTION 

C I LAS- 1 3  (814) 

C I LAS- 1 3  (815) 

CHARGE  WT.  (gm) 

1.7004  PbN6 

1.7099  PbN6 

CHARGE  RADIUS  (cm) 

.508 

.508 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

KAFB  Alluvium 

KAFB  Alluvium 

TEST  BED  DENSITY  (gm/cc) 

1.622 

1.581 

MOISTURE  CONTENT 

nom.  4.5% 

nom.  4.5?' 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CRATER  DIAMETER  (cm) 

12.94 

13.8 

CRATER  RADIUS  (cm) 

6.47 

6.90 

MAX.  CRATER  DEPTH  (cm) 

1.85 

1.95 

CRATER  ASPECT  RATIO  (r/h) 

3.50 

3.54 

CRATER  VOLUME  (cc) 

58.7 

72.4 

LIP  DIAMETER  (cm) 

23.0 

20.0 

LIP  HEIGHT  (cm) 

.163 

.131 

LIP  VOLUME  (cm3) 

104. 

61.7 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS  FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

1.850 

1.950 

1.0 

1.588 

1.713 

2.0 

0.744 

0.988 

3.0 

0.619 

0.806 

4.0 

0.613 

0.631 

5.0 

0.269 

0.369 

6.0 

0.006 

0.081 

7.0 

-0.044 

0.006 

8.0 

-0.031 

-0.031 

9.0 

-0.088 

-0.125 

10.0 

-0.156 

-0.131 

11.0 

-0.163 

-0.119 

12.0 

-0.163 

-0.088 

13.0 

-0.119 

-0.075 

14.0 

-0.113 

-0.063 

15.0 

-0.081 

-C. 005 

16.0 

-0.056 

-0.031 

17.0 

-0.056 

-0.025 

18.0 

-0.044 

-0.025 

19.0 

-0.038 

-0.019 

20.0 

-0.031 

-0.013 

21.0 

-0.019 

0.000 

22.0 

-0.019 

0.000 

23.0 

-0.019 

-- 

24.0 

-0.006 

-- 

129 


Table  B5  (cont.)  KAFB  alluvium  1G  control  shot,  1.7  grams  lead  azide. 


SHOT  NUMBER 

UDRI-648 

DATE 

18  Nov.  77 

PURPOSE 

1G  Control 

CHARGE  DESCRIPTION 

CILAS-13  (816) 

CHARGE  WT.  (gm) 

1.7004  PbN6 

CHARGE  RADIUS  (cm) 

.508 

CHARGE  CONFIGURATION 

half  buried  sphere 

TEST  BED  MATERIAL 

KAFB  Alluvium 

TEST  BED  DENSITY  (gm/cc) 

1.600 

MOISTURE  CONTENT 

nom.  4.5‘" 

TEST  BED  GEOMETRY 

homogeneous 

CRATER  DIAMETER  (cm) 

11.68 

CRATER  RADIUS  (cm) 

5.84 

MAX.  CRATER  DEPTH  (cm) 

1.75 

CRATER  ASPECT  RATIO  (r/h) 

3.34 

CRATER  VOLUME  (cc) 

45.0 

LIP  DIAMETER  (cm) 

20.0 

LIP  HEIGHT  (cin) 

.144 

LIP  VOLUME  (cm3) 

71.8 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS 

RANGE  (cm) 

DEPTH  (cm) 

0.0 

1.750 

1.0 

1.450 

2.0 

0.713 

3.0 

0.581 

4.0 

0.425 

5.0 

0.113 

6.0 

-0.075 

7.0 

-0.006 

8.0 

-0.044 

9.0 

-0.081 

10.0 

-0.144 

11.0 

-0.144 

12.0 

-0.119 

13.0 

-0.113 

14.0 

-0.113 

15.0 

-0.075 

16.0 

-0.038 

17.0 

-0.025 

18.0 

-0.013 

19.0 

0.000 

20.0 

0.000 
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Figure  B5 .  Comparison  of  craters  646,  647  and  648. 
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Table  B6.  KAFB 

alluvium  1G  control  shots,  1 

265  grams  PETN. 

, 

H 

H 

SHOT  NUMBER 

UDRI-649 

UDRI-652 

DATE 

22  Nov.  77 

9  Dec.  77 

PURPOSE 

1G  Control 

1G  Control 

. 

CHARGE  DESCRIPTION 

CICS-1 . 265  (B-6) 

CICS-1. 265  (B-7) 

i 

CHARGE  WT.  (gm) 

.125  AgN3/1.2241  PETN 

.125  AgN3/1.2259  PETN 

* 

s 

CHARGE  RADIUS  (cm) 

.565 

.565 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

1 

TEST  BED  MATERIAL 

KAFB  Alluvium 

KAFB  Alluvium 

\ 

TEST  BED  DENSITY  (gm/cc) 

1.571 

1.587 

MOISTURE  CONTENT 

nom.  4.5% 

nom .  4.5% 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CRATER  DIAMETER  (cm) 

13.75 

12.88 

CRATER  RADIUS  (cm) 

6.88 

6.44 

MAX.  CRATER  DEPTH  (cm) 

2.40 

2.00 

CRATER  ASPECT  RATIO  (r/h) 

2.87 

3.22 

CRATER  VOLUME  (cc) 

103. 

66.6 

LIP  DIAMETER  (cm) 

16.0 

14.0 

LIP  HEIGHT  (cm) 

.131 

.075 

LIP  VOLUME  (cm3) 

65.8 

60.8 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS  FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

2.400 

2.000 

1.0 

2.244 

1.738 

2.0 

1.513 

0.844 

' 

3.  C 

0.956 

0.688 

4.0 

0.838 

0.606 

5.0 

0.613 

0.369 

6.0 

0.250 

0.081 

7.0 

-0.044 

-0.075 

8.0 

-0.131 

-0.050 

9.0 

-0.038 

-0.013 

10.0 

0.000 

0.013 

11.0 

-0.069 

0.019 

12.0 

-0.113 

0.006 

l 

13.0 

-0.056 

0.000 

1 

14.0 

-0.050 

-0.025 

fl 

15.0 

-0.069 

-0.050 

16.0 

-0.019 

-0.044 

17.0 

-0.019 

-0.031 

18.0 

-0.038 

-0.056 

19.0 

-0.025 

-0.056 

20.0 

-0.025 

-0.044 

Mf  HI  lllf  li  |  ||ii"J  “ 

■MBMt 

21.0 

-0.025 

-0.056 

BMlBi 

22.0 

-0.019 

-0.025 

PPiMIKl»r 

23.0 

-0.013 

-0.025 

24.0 

-0.019 

-0.025 

L _ ! _ 
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Figure  B6.  Comparison  of  crater  649  with  crater  652. 
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Table  87.  KAFB  alluvium  1G  control  shots,  4  grams  PETN 


SHOT  NUMBER 

UDRI-650 

UDRI-651 

DATE 

30  Nov.  77 

5  Dec.  77 

PURPOSE 

1G  Control 

1G  Control 

CHARGE  DESCRIPTION 

CICS-4  (BL-8) 

CICS-4  (BL-9) 

CHARGE  WT.  (gm) 

.125  AgN3/3.9509  PETN  . 

125  AgN3/3.9620  PETN 

CHARGE  RADIUS  (cm) 

.826 

.826 

CHARGE  CONFIGURATION 

half  buried  sphere 

half  buried  sphere 

TEST  BED  MATERIAL 

KAFB  A1 luvium 

KAFB  Alluvium 

TEST  BED  DENSITY  (gm/cc) 

1.555 

1.584 

MOISTURE  CONTENT 

nom .  4% 

nom.  4'' 

TEST  BED  GEOMETRY 

homogeneous 

homogeneous 

CRATER  DIAMETER  (cm) 

16.56 

17.0 

CRATER  RADIUS  (cm) 

8.28 

8.52 

MAX.  CRATER  DEPTH  (cm) 

2.86 

2.80 

CRATER  ASPECT  RATIO  (r/h) 

2.90 

3.04 

CRATER  VOLUME  (cc) 

196. 

175. 

LIP  DIAMETER  (cm) 

20.0 

20.0 

LIP  HEIGHT  (cm) 

.30 

.294 

LIP  VOLUME  (cm3) 

354. 

403. 

MEAN  CRATER  PROFILE 

AVERAGE  OF  8  RADIALS 

FOR  EACH  CRATER 

RANGE  (cm) 

DEPTH  (cm) 

DEPTH  (cm) 

0.0 

2.850 

2.800 

1.0 

2.863 

2.663 

2.0 

2.425 

2.156 

3.0 

1.719 

1.400 

4.0 

1 .213 

1.038 

5.0 

1.038 

0.844 

6.0 

0.763 

0.650 

7.0 

0.388 

0.413 

8.0 

0.075 

0.156 

9.0 

-0.200 

-0.131 

10.0 

-0.300 

-0.294 

11.0 

-0.250 

-0.288 

12.0 

-0.163 

-0.219 

13.0 

-0.119 

-0.125 

14.0 

-0.094 

-0.100 

15.0 

-0.056 

-0.100 

16.0 

-0.094 

-0.100 

17.0 

-0.100 

-0.119 

18.0 

-0.088 

-0.125 

19.0 

-0.081 

-0.106 

20.0 

-0.113 

-0.088 

21.0 

-0.075 

-0.075 

22.0 

-0.075 

-0.075 

23.0 

-0.081 

-0.094 

24.0 

-0.063 

-0.075 
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MATERIAL  PROPERTIES  OF  TEST  SOILS 


REPORT  OP 
LABORATORY  TESTS 
ON  SOIL  MATERIALS 
FOR 

BOEING  AEROSPACE  COMPANY 


DECEMBER  1977 


SHANNON  &  WILSON,  INC. 
GEOTECHNICAL  CONSULTANTS 
1105  NORTH  38TH  STREET 
SEATTLE,  WASHINGTON  98103 


LriECEDlNG  PAOS 


gi  t  UK -NOT  F1LWED 


137 


TABLE  OP  CONTENTS 


Page 

L 

PURPOSE  AND  SCOPE 

139 

n. 

TEST  PROCEDURES 

139 

m. 

TEST  RESULTS 

140 

A.  Presentation 

140 

B.  Interpretation 

140 

Table 

Cl 

LIST  OF  TABLES  AND  FIGURES 

Summary  of  Materials  and  Test  Specifications 

141 

Figure 

Cl 

Modeling  Clay  Test  Results 

142 

Figure 

C2 

Sawing  Sand  Test  Results 

143 

Figure 

C3 

Flintshot  Sand  Test  Results 

144 

Figure 

C4 

Banding  Sand  Test  Results 

145 

Figure 

C5 

KAFB  Alluvium  #2  Test  Results 

146 

Figure 

Cb 

KAFB  Alluvium  #4  Test  Results 

147 

138 


REPORT  OP 
LABORATORY  TESTS 
PERFORMED  OH 
SOIL  MATERIALS 

i.  run  POSE  AND  SCOPE 


This  report  presents  the  results  of  a  series  of  unconsolidated-undrained  (UU), 
triaxial  compression  tests  performed  on  each  of  six  soil  materials  provided  by  Boeing 
Aerospace  Company.  Three  UU  tests  were  performed  on  each  material  type  to 
determine  the  shear  strength  parameters  of  the  material.  The  samples  were  designated 
as  follows:  KAFB  Alluvium  #2,  Banding  Sand,  Sawing  Sand,  KAFB  Alluvium  #4,  Molding 
Clay  and  Flintshot  Sand. 


II.  TEST  PROCEDURES 

Three  test  specimens  were  prepared  from  each  of  the  six  material  types.  Test 
specimens  of  the  five  .-.and  and  alluvium  material  types  were  prepared  by  compacting 
the  material  to  a  density  specified  by  Boeingfsee  Table  Cl  for  a  summary  of  specified 
densities).  Compaction  to  the  specified  density  was  achieved  by  vibrating  a  known 
weight  of  material  into  a  cylindrical  mold  of  known  volume.  The  mold  was  lined  with  a 
thin  rubber  membrane  that  remained  in  place  around  the  compacted  test  specimen 
throughout  the  test  in  order  to  protect  the  material  from  the  water  used  as  a  confining 
m-dium  during  the  test.  Test  specimens  from  the  "modeling  clay"  material  type  were 
prepared  by  trimming  the  bulk  sample  into  three  cylindrical  specimens  and  then 
encasing  each  specimen  in  a  thin  rubber  membrane  to  protect  them  from  the  confining 
medium. 

After  each  test  specimen  was  prepared,  it  was  mounted  in  a  triaxial  test 
chamber  and  subjected  to  a  specified  triaxial  comfin:ng  pressure,  then  it  was 
immediately  sheared  under  strain-controlled  conditions,  while  maintaining  a  constant 
confining  pressure  and  without  allowing  drainage.  For  each  material  type,  triaxial 
confining  pressures  of  100,  200  and  400  psi  were  used  on  the  first,  second  and  third  test 
specimens,  respectively,  as  specified.  Water  was  used  as  the  confining  medium  for  all 


tests. 


W-3306-01 


Simultaneous  readings  of  load  and  deformation  were  obtained  at  regular  time 
intervals  throughout  the  shearing  period.  Deviator  stress  and  axial  strain  values  were 
then  computed  from  these  readings. 

III.  TFSTR'^HI.TS 


A .  Presi  21 1  a  Upn 

The  results  of  each  of  the  six  series  of  UU  tests  are  presented  in  Figures 
Cl  thru  C6.  Included  in  each  figure  are  plots  of  deviator  stresses  versus  axial  strains, 
Mohr  circles,  pertinent  specimen  dimensions  and  test  parameters  for  each  of  the  three 
tests  in  the  series.  In  addition,  sketches  of  each  specimen's  failure  mode,  failure 
criteria,  and  sample  classification  are  included. 

B.  Interpretation 

The  results  of  these  triaxial  compression  tests  may  be  used  to  estimate  the 
parameters  of  cohesion,  c,  and  angle  of  internal  friction,  t  ,  in  the  Mohr-Conlomb 
equation  for  shear  strength  under  unconsolidated,  undrained  conditions.  This  equation 
relates  shear  strength  to  triaxial  confining  pressure  by: 

t=  c  +  N  tan  i> 

where:  t  =  shear  strength 

and  N  =  principal  stress 

The  parameters,  c  and  t ,  are  obtained  from  the  line,  drawn  tangent  to  the  three 

Mohr  circles  of  a  test  series,  which  represents  the  Mohr-Conlomb  equation.  The  Mohr 

circles  are  drawn  with  centers  at  +  ^°1  ~  q3^  on  the  principal  stress  axis  and  with 

radius  ,°1  °3*  ;  where  oj  is  the  triaxial  confining  pressure  and  (o,  -o„)  is  the  deviator 

\ — 2 - /  1  J 

stress.  Tne  values  of  deviator  stress  used  to  construct  the  Mohr  circles  are  obtained 

from  the  plots  of  deviator  stress  versus  axial  strain  based  on  the  failure  criteria  used 

for  the  test  series. 
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SAMPLE 


Modeling  Clay 
Modeling  Clay 
Modeling  Clay 

Sawing  Sand 
Sawing  Sand 
Sawing  Sand 

Flintshot  Sand 
Flintshot  Sand 
Flintshot  Sand 

Banding  Sand 
Banding  Sand 
Banding  Sand 

KAFB  Alluvium 
KAFB  Alluvium 
KAFB  Alluvium 

KAFB  Alluvium 
KAFB  Alluvium 
KAFB  Alluvium 


TABLE  Cl 

SUMMARY  OF 
MATERIALS  AND 
TEST  SPECIFICATIONS 


SPECIFIED  DENSITY 
3 

gms/cc _ lbs/ft 


SPECIFIED 

CONFINING 

PRESSURE  LABORATORY 

psi  TEST  NUMBER 


— 

— 

1.78 

110.1 

1.78 

110.1 

1.78 

110.1 

1.78 

110.1 

1.78 

110.1 

1.78 

110.1 

1.68 

104.8 

1.68 

104.8 

1.68 

104.8 

*2 

1.61 

100.5 

#2 

1.61 

100.5 

#2 

1.61 

100.5 

*4 

1.61 

100.5 

*4 

1.61 

100.5 

#4 

1.61 

100.5 

100 

UU-501 

200 

UU-502 

400 

UU-503 

100 

UU-301 

200 

UU-302 

400 

UU-303 

100 

UU-601 

200 

UU-602 

400 

UU-603 

100 

UU-201 

200 

UU-202 

400 

UU-203 

100 

UU-101 

200 

UU-102 

400 

UU-103 

100 

UU-401 

200 

UU-402 

400 

UU-403 
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TRIAXIAl  COMPRESSION  TEST 
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20.0 


AO  .  0 


120.0  lAO.O 
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60 . 0  30.0  100.0 


Te  s  t  No . 


Con  fining 

Pressure,  tsf 
H  t .  ,  In. 

D  i  a  .  ,  in. 

Ht./Dia.  Ratio 
Dry  Unit  Wt.,  pcf 
Water  Content,  I 
Before  Test 
After  Test 


*g  .  t  S  t  ra i n /M i n  . 
jmpac t i on  ,  X  I 
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1  A  .  VJ 

28.73 

6.09 

6. 3'J 

6.65 

3.0  1 

3.00 

3.00 

2.23 

2.29 

2.23 

92.0 

93.  A 

92  .  3 

.  U 

A  .  A 

A.  A 

A  .  A 

A  .  A 

A  .  A 

.706 

.  721 

.699 

CLASSIFICATION: 

Light  brown,  sandy  SILT. 


FIGURE  C6 
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APPENDIX  D 

STEREO-CAMERA  DYNAMIC  PHOTOGRAMMETRY 


Analytical  photogrammetry  can  be  applied  to  obtain  spatial  measurements 
of  explosive  crater  formation  at  the  time  of  formation  while  under  large 
centripetal  acceleration  in  an  operating  centrifuge.  Close-range  photogrammetry 
techniques  have  been  developed  for  small  craters  in  granular  soils  using 
nonmetric  cameras.  The  surface  area  and  the  volume  of  a  calibrated  object  space 
control  frame  photographed  under  static  1-G  conditions  have  beer  determined.  In 
addition,  a  contour  map  was  prepared  for  a  typical  crater  also  photographed 
under  static  conditions.  These  show  a  high  degree  of  accuracy.  In  preparation 
for  dynamic  tests,  cameras  have  been  made  to  operate  under  the  high  centripetal 
accelerations  of  the  rotor  hub  during  centrifuge  operation.  It  remains  to  be 
shown  what  accuracy  degradation,  if  any,  results  from  these  centrifugal  loads 
acting  on  the  camera  system. 

DESIGN  CRITERIA 

The  conventional  profilometer  method  of  post-test  measurement  can  be 
employed  only  for  stable  craters.  For  fluid  craters  the  topological  mapping 
must  be  obtained  dynamically  under  high  centripetal  accelerations.  It  was  the 
purpose  of  this  study  to  develop  a  relatively  rapid  and  accurate  method  for 
obtaining  such  data  through  the  application  of  photogrammetric  techniques. 

Due  to  the  fact  that  clearance  volume  within  the  arm  of  the  centrifuge 
is  limited  to  a  length  of  available  140  cm  with  a  38-cm  width  and  a  15-cm  depth, 
the  maximum  distance  between  the  camera  film  plane  and  the  specimen  is  100  to 
140  cm  depending  upon  the  size  of  the  test  sample.  The  maximum  base  of  the 
camera  station  is  ?5  cm  and  the  height  of  the  cameras  can  not  exceed  15  cm. 
Since  the  photographs  are  to  be  taken  of  specimens  undergoing  centripetal 
accelerations  up  to  600  G,  the  vibration  of  the  camera  system  must  also  be  taken 
into  consideration  in  order  to  obtain  quality  images. 

PHOTOGRAMMETRIC  SYSTEMS 


Since  there  are  no  terrestrial  cameras  with  fixed  base  that  meet  these 
criteria,  two  motor-driven  Nikon  F2  35-rm  cameras  were  considered  in  conjunction 
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with  analytical  concepts.  A  control  frame  was  built  around  the  specimen  to 
provide  an  image  on  the  photographs  which  can  be  measured,  allowing  the  camera 
positions  and  orientation  relative  to  the  frame  to  be  computed.  This  allows 
calibration  for  measuring  three  dimensional  coordinates  of  the  crater. 

The  two  cameras,  each  having  a  focal  length  of  55  mm,  are  mounted  on  a 
fixed  camera  base  of  21  cm  within  the  arm  of  the  centrifuge  as  shown  in  Figures 
D1  and  D2.  The  fixed  control  frame  with  20  reference  points  at  five  different 
levels  machined  into  the  soil-specimen  container  is  shown  in  Figure  D3.  The 
target  diameter  of  the  object-space  control  points  was  selected  using  eq.  D1 . 

T  _  Dm  (D1 ) 

T  =  j- 

where 

T  =  target  diameter  on  control  frame 

D  =  distance  between  camera  and  specimen 

f  =  focal  length 

m  =  diameter  of  micrometer  measuri ng-dot  on  comparator 

In  order  to  fit  the  various  measuring  instruments,  two  different  targel 
diameters,  0.05  cm  and  0.10  cm,  were  drilled  into  alternating  levels  of  the 
control  frame  which  produce  film  image  sizes  of  approximately  25  p  and  50  u , 
respecti vely. 

The  on-board  location  of  the  cameras  was  chosen  to  avoid  the  complex 
optical  geometry  which  would  have  been  necessary  with  a  stationary  mount.  In 
addition,  the  large  tangential  velocity  of  the  rotor  tip  would  have  required  a 
strobe  light  with  a  sub-microsecond  duration. 

The  cameras  were  mounted  as  close  to  the  rotor  axis  as  the  six-inch- 
diameter  hub  would  allow.  The  actual  location  of  the  cameras  produced  an 
average  acceleration  in  excess  of  60  G  on  each  camera  at  the  maximum  centrifuge 
speed  of  620  rpm.  Tests  were  performed  photographing  high  contrast  resolution 
charts  mounted  on  the  rotor  tip.  The  cameras  performed  sat i sfactor i ly  with  no 
appreciable  loss  of  resolution  up  to  maximum  centrifuge  speed.  The  shutter 
mechanism  worked  properly  on  both  cameras  throughout  the  entire  speed  range. 
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D3.  Control  frame  and  control  points. 


The  cameras  were  operated  in  the  semiautomatic  mode  using  the  motors  to 
advance  the  film.  At  600  rpm  the  motor-drive  mechanism  of  the  right-hand  camera 
started  to  advance  the  film  somewhat  erratically.  The  location  at  a  radius  of 
approximately  19  cm  corresponds  to  77  G.  This  can  be  remedied  by  mounting  the 
right-hand  camera  upside  down  bringing  the  wind  mechanism  more  inboard.  The 
camera  motors  were  powered  through  the  centrifuge  slip  rings  using  a  standard 
laboratory  16-Vdc  supply. 

In  the  automatic  mode  at  5  frames  per  second,  the  shutter  speed  (set  at 
1  /500  of  a  second)  went  out  of  calibration  at  approximately  300  rpm  corres¬ 
ponding  to  15  to  20  G  acting  on  the  cameras.  No  further  centrifuge  tests  were 
performed.  A  slower  shutter  speed  with  an  electronic  flash  may  provide 
satisfactory  operation  at  higher  rpm  in  the  automatic  framing  mode.  In  general, 
the  operation  of  the  cameras  exceeded  expectations  and  it  is  felt  that  they  will 
provide  satisfactory  photographs  for  the  intended  dynamic  photogrammet ry 
application.  Figure  D4  shows  typical  dynamic  photos. 

METHODOLOGY 

Analytical  photogrammetry  is  not  a  new  concept,  having  been  thoroughly 
developed  by  mathematicians  prior  to  the  1 930 1 s .  However,  the  laborious 
computations  inspired  development  of  analogical  instruments  which  make  a  simple 
graphic  solution.  The  speed  with  which  present  day  computers  perform  computa¬ 
tions  has  revitalized  the  mathematical  analytical  methods  which  permit 
correction  of  systematic  errors  not  correctable  with  the  mechanical  analogical 
system.  ^  ^ 


Analytical  photogrammetry  uses  data  measured  stereoscopical ly  using  a 
comparator  in  a  coordinate  system  with  the  origin  at  the  principal  point  based 
on  camera  fiducial  marks.  These  data  are  transformed  to  photo  coordinates  using 
a  linear  transformat  ion  to  change  the  scale  and  rotation. 


11  R  +  ai?  H  +  xQ 

(02. 1) 

n  h  -  ai?  R  +  y0 

(D2.2) 
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is  of  craters  under  the  acceleration 
id  1 00G  (below). 


The  transformed  coordinates  of  a  point  P  in  the  photo  coordinates  system  are  x 
and  y.  R  and  H  are  the  coordinates  of  point  P  expressed  in  the  comparator 
system,  a-^  and  a-|2  are  unknown  coefficients  of  the  transformation,  and  x0  and 
yQ  are  the  unknown  translation  terms  in  x  and  y. 

Since  the  Nikon  camera  is  not  equipped  with  fiducial  marks,  comparator 
coordinates  of  the  four  corner  points  were  measured  and  used  in  eq.  0?,  with  a 
least-squares  method,  to  determine  the  transformation  coefficients.  These  were 
then  used  on  all  the  points  in  each  photo  for  transforming  to  the  photo 
coordinates  system.  The  camera  station  for  each  photo  was  then  determined  using 
the  coordinates  of  the  object-space  control  points  and  their  corresponding  photo 
coordinates,  using  collinearity  equations^’^  as  follows: 


(E 

tr 


Eo)m 


n 


(N  -  Nq)  m12  +  (Z  -  ZQ)  m13 


m31  +  V  m32  +  V  m33 


■[ 


(E  -  Eq)  m2]  +  (N  -  Nc)  m 


(E  -  E„)  m31  +  (» 


22 


+  (Z  -  Z0)  m 


N0)  m32  +  (Z  -  Z0)  m 


] 

33  "I 

33  J 


(D3.1) 


(03.?) 


where  x  and  y  are  the  photo  coordinates  of  the  control  point  whose  ground 
coordinates  are  E,  N,  and  Z.  EQ,  N  ,  ZQ  are  the  coordinates  of  the  camera 
station,  f  is  the  focal  length  of  the  camera ,  •  are  elements  of  the  rotational 

orthogonal  matrix  consisting  of  direction  cosines  or  of  the  exterior  orientation 
elements  (<t>,  w ,  k)  of  the  presently  unknown  camera  station.  The  above  equations 
can  be  written  in  the  form  of  the  collinearity  condition  using  determinant 
notation 


x 


I^E* 

x 

m2e* 


=  0  =  Fx 

=  0  =  Fy 


(IVJ.l  ) 


(r>4.?) 
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where 


and 


Mi  "  mil »  mi2’  mi3 


E*  = 


E  -  E, 


N  -  N 
Z  -  Z . 


o 


mll 

= 

cos  <|> 

cos  < 

m12 

= 

COS  4> 

sin  ic 

m1 3 

= 

sin  <t> 

m21 

= 

-si  n  w 

sin  <t> 

COS  PC 

-  COS  ID 

sin  ic 

m22 

= 

-sin  a) 

sin  $ 

sin  ic 

+  COS  CD 

COS  PC 

m23 

= 

cos  $ 

si n  o) 

m31 

= 

-cos  u 

sin  if 

COS  PC  + 

sin  id  cos 

PC 

m32 

= 

cos  w 

sin 

sin  ic  - 

sin  cd  cos 

IC 

m33 

= 

COS  co 

COS  if 

Si  nee 

the  number  of  control  points 

exceeds 

least-squares  adjustment  is  used.  The  observation  equations  are: 


Vx  =  Fx  +  dFx 


V  =  F  +  dF 

y  y  y 


(06.1) 

(06.2) 


where  F's  are  composed  of  initially  estimated  values  and 


3  E  3  F 

dFx  =  ir  dEo +  w- 
o  o 


dN. 


3  F 

+  TT  dZo 

0 


3  F 


3id 


3  F 
_ ) 

3 


d$ 


3  r 


3  ic 


(06) 


The  operator  d  denotes  corrections  to  the  initial  approximat  ions ,  and  the  V’s 
are  residuals  of  the  photographic  measurement.  Fq.  06  is  substituted  into  eq. 
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D5 ,  which,  after  1  ineari zation  using  a  Taylor  expansion  neglecting  all  second 
and  higher  order  terms,  becomes 


Vx  =  b]dE0  +  b2dN0  +  b3df  +  b4d^  +  bjd*  +  bgd<  +  lx 

where 

lx  =  Fx  -  x 


bl = 

-f 

FLE* 

M]  E* 

(m3e*)2 

J 

a(M3E*) 

3t  0 

a(M.E*) 

h  = 

-f 

¥* 

MjE* 

D2 

(M3E*)2 

3 (M3E*) 

a(M-|E*) 

3No 

3N 

0 

h  =: 

-f 

¥* 

M^E* 

b3 

(M3E*7 

3(M3E*) 

»(M1E*) 

37o 

h  = 

-f 

M3E* 

M^E* 

°4 

? 

(M  E*) 

3(M3E*) 

3(M1E*) 

3w 

9uj 

h  = 

-f 

M3E* 

M^E* 

b5 

(m3e*)2 

9(M3E*) 

3(M]E*) 

— w~ 

b6  = 

-f 

M^E* 

(«3e*)2 

3(M3E*) 

»(m7e*) 

- 

3«r 

(07 ) 
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A  similar  equation  is  formpd  for  the  dFy.  Equations  D5  are  solved  for  the 
differential  corrections,  which  are  then  added  to  the  approximat ions  to  updatp 
them.  The  iteration  is  continued  until  the  corrections  become  negligibly  small. 

The  positions  of  the  cameras,  as  well  as  the  orientation  elements,  are 
determined  through  these  resections.  The  location  of  the  points  observed  on  the 
crater  are  determined  by  space  intersections.  The  same  mathematical  concept  is 
then  repeated  with  the  only  difference  being  the  coordinates  of  camera  stations 
(Eq,  Nq,  Z  )  now  known  as  E,  N,  and  Z  space  coordinates.  More  detailed  working 
equations  can  be  found  in  References  D1  through  D4.  The  generalized  flow  chart 
for  the  above  computations  is  shown  in  Figure  05. 

CONTROL  FRAME  CALIBRATION  RESULTS 

Four  ground  control  points  were  used  for  the  determination  ot  the 
camera  stations  and  the  orientation  elements.  The  other  16  ground  control  points 
were  used  as  check  points  for  evaluation  of  the  accuracy.  The  results  of 
analysis  of  photographs  taken  under  static  1-G  conditions  of  the  checl  points 
are  shown  in  Table  D1 .  As  can  be  seen  from  this  table,  the  coordinates  of  the 
check  points,  from  the  photogrammet ric  method,  when  compared  to  the  results  of 
the  actual  measurements  differ  by  only  0.0?  cm  and  0.04  cm  in  the  two  horizontal 
directions  and  0.14  cm  in  the  vertical  direction. 

The  surface  area  and  volume  of  the  control  frame  (as  well  as  of  the 
crater)  can  be  obtained  by  using  the  E,  N  and  Z  coordinates,  which  are 
determined  by  the  space  intersection.  The  total  surface  area  may  be  computed 
from  the  following  equation:^ 


s .  =  J 

i  =  n 

E 

,/( F?  . 

J  ? 

i  =  1 

y  i  ,j 

. ,  N  .  . , 
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i  ,J 

The  total 
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Figure  D5. 


The  generalized  flow  chart  for  photogranrnetric  computation 


Table  D1 .  The  results  of  distance  determination  between  check 
points  on  the  control  frame. 


No. 

E  (x) 
cm 

N  (y) 
cm 

Z 

cm 

Gi  ven 

Comp. 

Diff. 

Given 

Comp 

Diff. 

Gi  ven 

Comp . 

Diff. 

11 

100.00 

99.98 

+0.02 

127.33 

127.27 

+0.06 

104.00 

103.82 

+0.18 

12 

100.00 

100.00 

0.00 

126.33 

126.23 

+0.10 

103.00 

102.76 

+0.24 

13 

100.00 

100.00 

0.00 

125.33 

125.26 

+0.07 

102.00 

101.77 

+0.22 

14 

100.00 

100.00 

0.00 

124.33 

124.28 

+0.05 

101.00 

100.87 

+  0.13 

*15 

100.00 

100.00 

0.00 

123.33 

123.31 

+0.02 

100.00 

99. °2 

+0.08 

21 

127.33 

127.27 

+0.06 

100.00 

100.00 

0.00 

104.00 

104.04 

-0.04 

22 

126.33 

126.31 

+0.06 

100.00 

100.01 

-0.01 

103.00 

103.12 

-0.12 

23 

125.33 

125.32 

+0.02 

100.00 

100.01 

-0.01 

102.00 

102.10 

-0.10 

24 

124.33 

124.34 

-0.01 

100.00 

100.12 

-0.02 

101.00 

101.16 

-0.16 

*25 

123.33 

123.36 

-0.03 

100.00 

100.00 

0.00 

100.00 

100.07 

-0.07 

31 

100.00 

100.00 

0.00 

72.67 

72.77 

-0.10 

104.00 

103.84 

+0.16 

32 

100.00 

100.01 

-0.01 

73.67 

73.71 

-0.04 

103.00 

102.76 

+0.24 

34 

100.00 

100.00 

0.00 

75.67 

75.71 

-0.04 

101.00 

100.86 

+0.1-1 

*35 

100.00 

100.00 

0.00 

76.67 

76.69 

-0.02 

100.00 

99.93 

+  0.07 

41 

72.67 

72.73 

-0.06 

100.00 

100.02 

-0.0? 

104.00 

103.91 

+0.09 

42 

72.6/ 

72.73 

-0.03 

100.00 

100.02 

-0.02 

103.00 

102.9° 

+0.09 

43 

74.67 

74.67 

0.00 

100.00 

100.01 

-0.01 

102.00 

102.03 

-0.03 

44 

75.67 

75.65 

+0.02 

100.00 

100.01 

-0.01 

101.00 

101 .10 

-0.10 

*45 

76.67 

76.65 

+0.02 

100.00 

100.00 

0.00 

100.00 

100.08 

-0.08 

Standard  error 

+0.02 

0.04 

0.1-3 

‘Used  as  control  points  for  resection. 
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1 


f  (E,N)  dEdN  = 


i=m 

s/sif  pi  (E,N)  dEdN 
i=l 


(P9) 


where  the  function  pi(E,N)  is  a  polynomial  of  fixed  degree  N,  i.e.,  Z  =  f  (E,N). 
With  this  approach,  the  surface  area  and  the  volume  of  the  control  frame  were 
determined  to  have  values  as  shown  in  Table  D2.  It  was  found  from  this 
experiment  that  a  0.23%  error  in  surface  area  was  obtained  and  a  0.62%  error  was 
observed  in  the  volume. 


Table  D2.  Results  of  the  surface  area  and  volume 
determination  for  the  control  frame. 


Sector 

Surface 

2 

cm 

Vol ume 

3 

cm 

Given 

Comp. 

Diff 

• 

Di  ff% 

Given 

Comp. 

Diff. 

Di  ff% 

1 

1089 

1088 

1 

0.09 

1089 

1074.4 

14.6 

1  .3 

4 

1387 

1382 

5 

0.36 

1387 

1379 

8.4 

0.56 

5 

1495 

1499 

4 

0.3 

1495 

1495 

0.0 

0.0 

Mean 

0.23% 

0.62% 

CRATER 

RESULTS 

Using 

the  control 

frame  calibration  of  the  photogrammetry  system,  a 

contour 

map  was 

prepared 

for 

a 

typical 

crater  photographed  statically.  This  is 

shown  in  Fig. 

06.  Depths 

on 

this  map  compare  favorably  with  those  of  the 

average 

crater 

contour 

shown 

in  Figure  D7,  which  was 

determined  using  eight 

radial  profilometer  measurements .  The  sterpophoto  pair  is  also  shown  in  Figure 
D7. 
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DEPTH  (cm) 


CONCLUSIONS 


Analytical  photogrammetry  has  a  great  potential  for  the  precise  mea¬ 
surements  of  explosively  formed  craters  under  centripetal  acceleration.  The 
present  experiment  gives  an  upper  limit  on  accuracy  based  on  a  static  calibra¬ 
tion.  The  surface  and  the  volume  of  the  control  frame  were  obtained  by  photo- 
gramnetry  with  accuracies  of  0.03  cm  in  the  horizontal  direction,  0.14  cm  in 
the  vertical  direction,  0.237-  in  surface  area  and  0.627  in  volume. 

These  results  were  obtained  from  the  photographs  taken  with  a 
nonmetric  camera  producing  no  fiducial  marks  on  the  photo  and  without 
correcting  for  lens  distortion.  It  is  expected  that  these  results  could  be 
considerably  improved  with  the  installation  of  fiducial  marks  on  the  present 
cameras.  A  calibration  procedure  should  be  used  to  determine  thp  true  focal 
length  of  the  camera  and  to  correct  for  lens  distortion  as  well  as  film 
distortion  under  high  gravity  forces. 

The  stability  of  the  system  under  gravity  and  with  continued  use  has 
not  been  conclusively  tested.  However,  from  all  indications  it  would  appear 
that  variations  in  centripetal  accel erations  would  have  only  a  small  effect  if 
exposure  as  well  as  processing  is  done  under  controlled  circumstances.  The 
residuals  in  the  analysis  of  the  check  control  points  suggest  that  an  error- 
existed  in  the  prior  measurements  of  the  control  frame  used  as  the  given 
reference  dimensions  or  the  establishment  of  thp  control  frame  may  contain  some 
systematic  errors,  both  of  which  can  be  corrected. 
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ATTN:  R.  Plebuch 
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ATTN:  E.  Wong 
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